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ABSTRACT 

The  purpose  of  this  thesis  is  the  investigation  of 
mathematical  methods  for  treating  collisions  of  light  nuclei. 

Our  understanding  of  nuclear  forces  is  derived  largely 
from  the  study  of  simple  nuclear  systems.   Bound  nuclei 
consisting  of  two  to  four  particles  have  been  theoretically 
treated  with  considerable  accuracy  and  yielded  many  valuable 
results.   Another  important  source  of  information  are  nuclear 
collisions,   Hov/ever,  up  to  the  present,  only  two  particle 
collisions  have  been  theoretically  treated  in  a  systematic 
and  reliable  manner. 

This  thesis  is  concerned  with  collisions  involving  any 
number  of  nucleons.   The  problem  is  to  calculate  the  scattering 
properties  of  nuclei  from  assumed  nuclear  interactions.   These 
theoretical  results  can  then  be  checked  against  experimental 
data, thus  providing  a  test  for  the  validity  of  the  interactions 
which  were  postulated.   As  the  wave  function  of  a  nuclear 
n  -  body  system  extends  over  5n  dimensions  and  is  highly 
irregular,  on  account  of  the  great  strength  of  nuclear  forces, 
one  cannot  hope  to  determine  it  with  great  exactitude  for 
systems  involving  more  than  two  particles •   For  this  reason 
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we  have  developed  variational  methods  v/hich  yield  the 
scattering  properties  of  nuclei  with  considerable  accuracy 
without  requiring  an  exact  knowledge  of  the  wave  function* 

The  wave  function  is  assumed  to  satisfy  a  Schroedinger 
equation  in  v/hich  nuclear  forces  are  described  by  two-particle, 
velocity-independent  potentials.   V/hen  the  colliding  nuclei 
are  far  apart,  the  form  of  the  wave  function  is  known:   It 
consists  of  a  product  of  functions  describing  the  internal 
states  of  the  nuclei,  and  a  free  particle   function,  P(r  ),  of 
the  relative  position  vector  r*  of  their  centers  of  gravity. 
To  predict  the  experimental  scattering  properties  of  the 
nuclei  the  detailed  characteristics  of  F( f  )   must  be 
determined  from  the  assumed  interactions. 

If  one  considers  collisions  in  which  the  nuclei  have 
a  definite  angular  momentum  of  relative  motion, Z  ,  the 
function  F(  r  )  has  the  asymptotic  form 

where  k  is  the  de  Broglie  wave  number  of  relative  motion, 
A  -s.    \y' \   $   S*  is.  the  angular  function  and  n^  the  phase 

shift  corresponding  to  the  angular  momentum  !L  •      The  phase 

shift   0^  determines  the  scattering  cross  section. 

Similarly,  the  function  P(r  )  corresponding  to  a 
plane  wave  of  one  nucleus  incident  on  another  has  the 

asymptotic  form 
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Here  1^,1   "••iki  ^  k,  k,  is  in  the  direction  of  incidence 
and  f(k^,  k)  is  the  scattering  amplitude  in  the  direction 
of  k. 

We  seek  to  determine  the  phase  shifts  and  scattering 
amplitudes,  given  the  interaction  between  the  nucleons.   In 
collisions  of  two  nucleons  the  wave  equation  separates  and 
its  solution  presents  no  essential  difficulty.   But  when  the 
colliding  nuclei  are  complex,  the  Schroedinger  equation  can 
no  longer  he  separated.   Hence  an  exact  solution  is  not 
possible  and  one  must  take  recourse  to  approximation  methods. 

We  have  followed  two  lines  of  attack,  both  of  which 
have  been  successful  in  bound  state  problems:   The  Rayleigh- 
Ritz  variational  method  of  undetermined  coefficients  and  the 
variation-iteration  method. 

Rayleigh-Ritz  Method* 

The  principle  of  one  of  the  procedures  adopted  may  be 
illustrated  by  the  one  dimensional  Schroedinger  equation. 
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which  is  to  be  solved  subject  to  the  boundary  condition  that 
u(0)  -s.  0.  The  interaction  potential  V(x)  may  be  supposed  to 
vanish  for  x  ^  a,  so  that 

We  now  form  the  expression 
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and  consider  its  first  variation,  restricting  our  trial 
functions,  ^^  >  ^7  the  conditions  u.(0)  =•  0  and 

After  an  Integration  by  parts,  we  find 

Since  for  the  correct  solution  u,  1^0,    it  follows  that 
tan  yi     determined  from  the  equation 

is  stationary,  relative  to  any  admissible  variation  of  u. 

When  this  equation  has  been  re-written  in  a  homogeneous 
form,  one  can  use  a  trial  function  of  the  type 

on  the  right  hand  side  and  determine  the  c  *s  from  the 
stationary  property  of  tan  vi  •   As  in  the  case  of  bound 
state  problems  this  leads  to  a  determinantal  compatibility 
equation.   In  contrast  to  secular  equations  this  is  a  linear 
equation  for  the  determination  of  tan  ri    ,  regardless  of  the 
number  of  parameters* 

This  and  similar  methods  are  applied  to  the  following 
cases: 


A.  Two  Particle  Collisions. 

1.  Phase  Shifts,   Stationary  expressions  are  developed 
for  the  phase  shifts  in  the  case  of  central  forces  and  for 
the  proper  phase  shifts  of  the  scattering  matrix  in  the  case 
of  tensor  forces.   Some  numerical  calculations  show  that  an 
accuracy  of  a  fraction  of  one  per  cent  is  attainable  without 
undue  labour. 

2.  Scattering  Amplitudes.   The  variational  method 
developed  permits  an  application  of  the  Rayleigh-Ritz  method; 
it  also  leads  naturally  to  the  Born  approximations  and 
Schv/ingei* '3  variation-iteration  formulation  (  J.S.  Schwinger, 
unpublished  lectures  ). 

3.  Level  Width.   A  stationary  expression  for  the  width 
of  a  resonance  level  is  derived. 

B.  Collisions  of  Composite  Nuclei. 

1.  Neutron  Deuteron  Scattering.   This  case  is  treated 
in  some  detail,  as  a  typical  example  of  composite  collisions. 
A  preliminary  calculation  is  carried  out  for  the  scattering 
cross  section  at  zero  energy  and  reasonable  agreement 

(  within  about  ten  per  cent  )  is  obtained  with  the  experimental 
value  of  3.2  barns. 

2.  General  Composite  Collisions.   A  few  features  of 
the  Rayleigh-Ritz  method  in  rather  general  collisions,  which 
may  involve  excitation  or  transmutation,  are  discussed. 
Variational  principles  for  the  elements  of  the  scattering 
matrix  are  derived. 
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where 


Variation  Iteration  Method. 

Instead  of  employing  trial  functions  with  undetermined 
coefficients  In  the  stationary  expressions  for  the  phase  shifts 
and  scattering  amplitudes  one  can  also  use  Iterated  functions. 

Thus  In  the  one  dimensional  problem  the  n-th  Iterate 
Is  defined  In  terms  of  the  (n-l)-st  by  the  equation 

0 

Professor  Schwlnger  (  J.  S.  Schwlnger,  unpublished  lectures  ) 
has  applied  this  method  very  successfully  In  collisions  of 
two  particles.   Since  the  iteration  improves  the  wave  function 
very  rapidly  this  method  converges  extremely  well.   Schwlnger »s 
method  is  generalized  from  twa  particle  collisions  to 
composite  collisions.   The  persistence  of  the  internal 
interactions  in  the  asymptotic  part  of  the  configuration 
space  causes  considerable  difficulty  in  the  construction  of 
a  suitable  Green's  function.   Of  several  methods  which  are 
investigated  only  one  (  involving  a  Green's  function  in  three 
dimensions  ),  is  found  at  all  feasible.   However  even  this 
procedure  is  so  complicated,  that  unless  a  more  suitable 
version  of  the  variation-iteration  method  can  be  found,  the 
Raylelgh-Rltz  method  appears  to  be  more  useful  for  the  treatment 
of  composite  nuclear  collisions.  ^  ,  yO  af^*<^*^Y^ 
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I  INTRODUCTION* 

A  great  deal  of  Information  about  the  nature  of  nuclear 
forces  has  been  derived  from  a  comparison  between  experimental 
and  theoretical  studies  of  simple  nuclear  systems.   Among  the 
most  important  data  so  compared  are  the  binding  energy  and  the 
magnetic  dipole  and  quadrupole  moments  of  the  deuteron  ;  the 
binding  energies  of  H  ,  He  and  He    ;  and  the  scattering 
cross  sections  of  neutrons  and  protons  by  protons   • 

However,  while  the  binding  energies  of  nuclear  3-  and  4- 
body  systems  have  been  theoretically  computed  and  have  led  to 
some  very  important  conclusions  about  nuclear  forces  (e.g»  the 
forces  operative  between  neutrons   ),  no  very  satisfactory 
scheme  for  mathematically  treating  collisions  which  involve 
more  than  two  elementary  particles  has  so  far  been  given.   In 
v/hat  follows  some  promising  methods  to  solve  such  problems 
will  be  investigated  and  preliminary  results  presented. 

Theoretical  calculations  of  collisions  of  light  nuclei 
are  worthy  of  interest  for  the  following  reasons: 

(1)   They  may  confirm  or  weaken  our  confidence  in  existing 
nuclear  theories;  in  particular,  they  may  throw  light  on  the 
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question  whether  nuclear  forces  are,  partly,  many  particle 
forces. 

(2)   They  may  help  to  decide  between  certain  ambiguities 
(  such  as  the  exact  space  and  spin  dependence  of  nuclear 
forces  ),  which  the  data  at  present  available  do  not  settle. 

From  the  very  fact  that  more  than  two  intimately  inter- 
acting particles  are  involved  it  is  clear  that  collisions  of 
composite  nuclei  present  a  very  complicated  theoretical  problem. 
V/hile  it  is  encouraging  that  the  corresponding  bound  state 
problems  have  been  solved  with  a  reasonable  amount  of  effort  to 
within  a  few  per  cent,  it  may  not  be  out  of  place  to  mention 
the  additional  complications  of  the  scattering  problems.   The 
asymptotic  form  of  the  wave  function,  instead  of  having  to 
vanish  at  infinity,  must  correctly  describe  the  colliding 
particles.   Above  certain  energy  threshholds  it  must  exhibit 
the  possibility  that  part  of  the  kinetic  energy  of  collision 
may  be  used  to  raise  one  or  both  of  the  colliding  particles  to 
an  excited  state;  or  that  an  exchange  of  nucleons  may  take 
place,  so  that  the  resulting  nuclei  differ  from  the  colliding 
ones;  (e.g.  T^^P^-^  W-e,  -tv\*   );  or,  finally,  that  a  triple 
or  multiple  disintegration  may  take  place  (e.g.  La"*  -4  P   — ^ 
(  intermediate  stages  )  -— ^  JlHe  -^v\       ).   Interactions 
with  the  radiation-  and  electron- neutrino  fields  will  not  be 
considered  In  this  paper.   We  shall  also  exclude  from  our   dis- 
cussion processes  in  which  multiple  disintegrations  are 
energetically  possible  since  in  that  case  the  asymptotic  form 
of  the  wave  function  is  very  complicated  and  it  was  considered 


advisable  to  concentrate  first  on  a  solution  of  the  simpler 
problems. 


Experimental  Data. 

In  the  last  twenty- five  years  a  great  many  experimental 
investigations  of  collisions  of  light  nuclei  have  been  carried 
out.   Many  of  those  performed  before  1937  are  discussed  In 
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Bethels  and  Livingston's  reviev/  articles  on  "Nuclear  Dymanlcs" 
A  few  other  typical  references  are  given  below     •   Scattering 
experiments  with  neutrons  as  of  1947  have  been  compiled  by 
H.  H.  Goldsmith  et  al,   .   Quite  recently,  collisions  at  very 
high  energies  have  been  investigated  by  means  of  the  Berkeley 
cyclotron   • 

It  is  our  general  aim  to  make  possible  a  theoretical 
Interpretation  of  these  experiments.   However,  this  paper  Is 
devoted  mainly  to  the  development  of  theoretical  methods  and 
only  very  few  numerical  calculations  have  so  far  been  carried 
out.   Hence  we  shall  not  be  concerned  here  with  a  detailed 
analysis  of  the  experimental  data. 
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Previous  Theoretical  Work» 

Theoretical  work  on  more  complicated  scattering  problems 
h^s  been  mostly  of  a  qualitative  and  phenomenological  nature. 
It  has  been  largely  based  on  the  resonance  formula  of  Breit 
and  V/igner  which  exhibits  the  scattering  process  as  the 
temporary  formation  of  a  compound  nucleus  with  subsequent 
emission  of  the  incident  or  some  other  particle.   In  terms  of 
this  theory  it  has  been  possible  to  understand  the  general 
features  of  nuclear  scattering,  but  no  detailed  information 
about  nuclear  forces  can  be  derived  from  it. 

V/heeler   and  Heisenberg  have  analyzed  the  general 
framev;ork  of  collision  processes  without  going  into  any  detail 
of  the  interaction.   Their  theory  serves  two  purposes:  It 
elucidates  the  general  features  of  scattering  phenomena, 
regardless  of  the  nature  of  the  interaction,  and  it  represents 
the  first  stage  of  any  detailed  calculation  of  composite 
collisions,  based  on  some  specific  nuclear  theory. 

Furthermore,  as  an  approximation  to  the  correct  wave 
function,  \'i'heeler  has  used  functions  of  so-called  group 
structure,  which  have  been  employed  in  all  succeeding 
calculations.   To  determine  the  best  possible  functions  of 
this  sort,  he  made  use  of  a  variational  principle.   The 
contents  of  sections  1.6.  and  1.7.  are  an  extension  of  this 
method. 
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More  directly  related  to  the  present  paper  are  the  works 
of  Schiff   ,  Ochlal   ,  Flttgge   ,  Motz  and  Schwinger   , 
Buckingham  and  Massey    ,  and  K6cker  who  have  calculated 
scattering  cross  sections  In  n-d  and  p-d  collisions  from 
approximate  solutions  of  the  Schroedlnger  equation. 

These  calculations  are  all  based  on  a  wave  function 
with  group  structure,  that  is  a  function  of  the  form 
(0.1)   Vj  .   £  ^;  C^.  -P^^ 

where  ^^:  ,  ^'\   describe  two  colliders  1  and  j  when 
outside  their  range  of  interaction,  and  ^ii    describes 
their  relative  motion.   The  sum  is  taken  over  all  pairs  of 
nuclei  whose  formation  is  compatible  with  the  given  energy 
and  over  all  permutations  of  identical  particles,  with  the 
signs  adjusted  in  accordance  with  the  exclusion  principle* 
The  functions  V^-  are  then  chosen  so  that  ju  satisfies  the 
Schroedlnger  equation  as  well  as  possible.   A  function  of  this 
type  has  the  correct  symmetry  properties  and  asymptotic 
behaviour;  but  in  view  of  the  strength  of  nuclear  forces  one 
cannot  expect  it  to  approximate  too  well  the  correct  v;ave  function 
in  that  region  of  the  configuration  space  where  the  collidors 
Interact,  and  it  is  just  this  region  which  determines  the 
scattering  cross  section. 
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Of  the  calculation  mentioned  above,  those  of  Buckingham 
and  Masaey   '  are  the  moat  extensive.   Their  method  is  in 
essence  the  following: 

The  centre  of  gravity  motion  of  the  two  neutrons  and 
proton,  participating  in  an  n-d  collision,  is  separated  out 
and  there  remains  the  following  equation  for  the  wave  function 
in  the  relative  coordinates: 

(0.2)    \r^  V  -t  ^^^  ^  ^^^v\  ^  -  (^w  +  t^]  ^ , 

Here  T  is  the  kinetic  energy  operator  for  the  relative  motion 
of  the  three  particles;  ^ <\    i3  the  interaction  potential 
between  particles  i  and  j,  supposed  independent  of  the  position 
of  the  third  particle;  ^i^  is  the  energy  of  relative  motion  of 
the  colliding  neutron  with  respect  to  the  deuteron  in  the  system 
in  which  the  centre  of  gravity  of  all  three  particles  is  at  rest; 
and  K^  is  the  binding  energy  of  the  deuteron. 

Let  us  introduce  relative  coordinates  to  which  we  shall 
have  frequent  reference  in  the  following  pages.   If  K^  ,  »>x. 
and  K^  are  the  absolute  position  vectors  of  the  two  neutrons 
and  the  proton,  then  P^  ,  9^  ,  9,   are  defined  by 

and  similarly  f^    ,     T^  ,    f^     are  defined  by 


(0.4)   ^   , 
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2. 
2^ 


Fig.  1 
In  terms  of  these  coordinates  T  has  the  form 

'-'  r-  -  ^  IH .  \  %:) 

or  the  two  other  forms  obtained  by  replacing  1  by  2  or  3» 
M  Is  the  mass  of  a  nucleon. 

The  wave  function  J^  Is  a  function  of  the  six  relative 
coordinates  represented  by  o     and   T   (  or  any  two  linear 
combinations  of  these  )  and  of  the  spins  of  the  three  particles* 

Equation  (0»2)  Is  to  be  solved  subject  to  the  condition 
that  when  either   K.   or   f^  become  large,   ju   describes  a 
deuteron  and  a  free  neutron  colliding  with  lt»   Thus 
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where  Q  f §)   is  the  ground-state  wave  function  of  the  deuteron, 

k'^  IkU  C'-^^en  /hm)'/'-  ,     A,x  \r.l^  .x^^\r,\     ,   f  i3  the 

scattering  amplitude  and  'Y^and'W^are  appropriate  spin 
functions  which  differ  according  to  whether  we  consider  quartet 
or  doublet  scattering. 

Although  we  are  not  familiar  with  a  rlgourous  proof  that 
a  slx-dlmenslonal  partial  differential  equation  of  elliptic 
type  has  a  solution  which  Is  finite  everywhere  and  satisfies 
such  asymptotic  "boundary  conditions,  one  Intuitively  expects 
this  to  be  the  case.   However  a  reasonably  rlgourous  solution 
Is  undoubtedly  very  complicated,   Buckingham  and  Kassey,  to 
simplify  their  problem  therefore  assume  that  the  exact  solution 
can  be  every^vhere  well  approximated  by  a  function  of  the  form 


(0 


.7)  5  ,  ^l^s)^l^H,    +  ^t?0  ^  Wl^^' 


This  Is  equivalent  to  allowing  only  for  the  distortion  of  the 
free  neutron  function  (  since  F  Is  as  yet  undetermined  ),  but 
disregarding  the  distortion  of  the  deuteron  by  the  neutron. 

To  find  which  function  F  makes  ^    the  best  possible 
solution,  they  substitute  (0,7)  into  (0.2)  and,  using  the  wave 
equation  for  the  deuteron,  they  find  that  F  should  satisfy,  in 
the  case  of  doublet  scattering, 

l^M<-+^H  -  i^-hK  ^  m\  ^tM  ^uo 

(0.8) 


Here  w,  h,  m  and  b  represent  the  strengths  of  the  Wlgner, 
Helsenberg,  Majorana  and  Bartlett  forces. 

For  any  fixed  value  of  ^.  ,  (0*8)  could  in  principle  be 
solved  for  ^l^)   giving'  ?iL^y'      ^^^  ^^  there  is  no  reason  why  P 
should  be  independent  of  P     it  is  impossible  to  satisfy  (0.8) 
by  a  single  function  f'(**)»   This  is  just  an  expression  of  the 
fact  that  a  function  of  the  type  (0,7)  cannot  solve  the  wave 
equation  exactly. 

To  obtain,  as  it  were,  an  average  solution  for  P, 
Buckingham  and  Massey  multiply  (0.8)  by  CP(§A  and  integrate 
over  the  P  -space.   This  leads  to  an  integro-differential 
equation  for  P,  of  the  form 

(0.9)  ^^"-"^"^  ^^'^  '  [^^-^)   lL(_v)^(r) 

where  U  and  Q  involve  the  known  interaction  potentials  and 
the  deuteron  function.  This  equation  is  finally  solved  by 
numerical  methods. 

Previous  to  Buckingham  and  Massey,  Motz  and  Schwinger 
have  calculated  the  scattering  of  neutrons  by  deuterons  in 
the  zero  energy  limit.   The  essential  features  of  their 
calculation  are  similar  to  those  of  the  British  workers. 
Thus,  they  too  consider  the  deuteron  as  rigid  and  are  led 
to  an  integro-differential  equation  for  the  function  P. 
To  solve  it  they  replace  it  by  an  integral  equation  according 
to  the  scheme : 
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(0.10) 


P(0  -  e 


After  carrying  out  angular  Integrations  analytically  the 
remaining  integral  equation  is  solved  approximately  by 
replacing  it  by  a  set  of  simultaneous  linear  equations. 

Finally,  Hticker  has  tried  to  write  V[rj  as  a  power 
series  in  .\  and  determine  the  coefficients  to  best  fit  the 
Schroedinger  equation.   However  his  results  depend,  by  a 
factor  of  about  1.5-2.0  on  the  point  where  the  pov/er-series 
is  taken  to  pass  over  into  the  asymptotic  form  of  F. 

The  other  calculations  do  not  contain  any  essentially 
different  methods. 

The  different  authors  using  various  kinds  of  (  central  ) 
interactions  obtained  widely  different  results.   For  the 
total  cross  section  of  n-d  scattering  at  0-energy  they  are  in 
units  of  barns  as  follows:  Schiff   :  6.11  ;  Ochiai   :  3.2  ; 
Motz  and  Schwinger    :  4.57,  6.91  j  Buckingham  and  Massey   : 
3.0,  3.9  (  extrapolated  )  ;  HQcker    :  1.88,  2.95  .   The 
experimental  value  is  near  3.3  '   . 

Other  two-  and  three-particle  systems  do  not  exhibit  such 
a  strong  dependence  on  the  type  of  interaction,  at  low  energies. 
If  taken  at  their  face  value,  then,  these  results  would  • 
provide  a  useful  method  for  distinguishing  between  various 
kinds  of  interaction,  all  of  which  are  compatible  with  other 

27.   L.  J.  Rainwater  et  al.,  Phys.  Rev.  73,  733  (1948). 
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low  energy  phenomena.      • 

However,  such  a  procedure  would  be  without  firm  foundation 
in  view  of  the  unreliability  of  the  group  structure  function, 
on  which  all  calculations  are  based.   This  emphasises  the 
need  for  a  systematic  method  of  solution,  which  is  not 
restricted  to  any  particular  kind  of  function* 


Summary  of  the  Present  Paper, 

,  It  is  not  the  primary  purpose  of  this  paper  to  add 
another  approximate  calculation  to  those  already  published, 
but  rather  to  investigate  the  possibility  of  a  systematic 
procedure  for  solving  the  Schroedinger  equations  corresponding 
to  composite  scattering  processes.   It  was  natural  to  look 
for  a  variational  method,  as  this  has  been  the  chief  tool 
for  solving  the  problems  of  the  bound  3-  and  4-particle 
systems    '    and  has  recently  been  applied  by  Schwinger 
to  2-particle  scattering  problems. 

There  are  two  types  of  practical  variation  procedures: 
One  is  the  Rayleigh-Ritz  method   of  undetermined  coefficients; 
the  other  the  variation  iteration  method  employing  Green's  functions 
We  have  tried  both  methods  in  the  many  particle  scattering 
problem,  and  we  strongly  incline  to   the  opinion  that  the 

28.  N.  Svartholm,  Thesis,  Lund,  Hakan  Ohlssons,  Boktryckeri  (1945). 

29.  Unpublished  lectures. 

30.  W.  Rltz,  J.  reine  angew.  Math.  135  (1909). 
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Rayleigh-Ritz  method  la  the  more  useful. 

Using  this  method  we  have  made  exploratory  calculations 
on  2-  and  3-partlcle  scattering  problems,  some  with  the 
inclusion  of  tensor  forces.   On  the  basis  of  these  we  believe 
that  this  procedure  can  indeed  be  used  for  a  systematic 
Investigation  of  collisions  of  light  nuclei. 
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I.   THE  RAYLEIGE-RITZ  METHOD  IN  COLLISION  PROBLEMS, 


I»l.   General  Remarks • 

In  his  original  paper  on  wave  mechanics,  Schroedlnger 
derived  the  wave  equation  from  a  variational  principle 
analogous  to  Hamilton's  Principle  in  classical  mechanics. 
The  variational  formulation  of  the  equations  of  motion  of 
wave-fields  has  been  very  fruitful  not  only  in  theoretical 

n 

developments  of  field  theory  but  also  in  approximate 
solutions  of  complicated  quantum  mechanical  problems   • 

The  method  most  extensively  used  in  practical  problems 
has  been  the  Rayleigh-Ritz  method  of  linear  trial  functions. 
However  its  application  has  been  limited  to  systems  in  bound 
states,  v/hose  wave  functions  vanish  rapidly  at  infinity. 

Recently,  Heisenberg   has  pointed  out  the  correspondence 
between  the  phase  shift  energy  relation  of  a  collision  problem 
and  the  energy  spectrum  obtained  if  the  participating  particles 
are  enclosed  in  a  large  box. 

Consider,  as  the  simplest  example,  a  particle  of  mass  M 
under  the  action  of  a  short  range  scattering  centre.   The  wave 
equation  for  the  radial  part  of  the  spherically  symmetric 
solution  is 

31.  E.  Schroedlnger,  Ann.  d.  Phys.  79,  361  (1926). 

32.  See,  for  example,.  G.  V/entzel,  "Einftihrung  in  die 
Quantentheorie  <i.  We llenf elder",  Wien,  Franz  Deuticke 
(1943). 

33.  E.  Hylleraas,  Zs.  f.  Phys.  54,  347  (1929). 
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(1 


where   U  y\\   is    the    scattering  potential   and  E   the    total   energy. 
Asymptotically   the   solution  becomes 


(1,2)        U.      >       ^'-^      ^JIE       A    -f 


T(ej\ 


where  the  phase  shift  ^[^]   >  corresponding  to  the  energy  E 
of  the  incident  particle  must  be  found  by  a  detailed 
solution  of  (1,1)« 

Now  let  us  ask  for  the  energy  levels  of  the  same 
particle  in  the  same  potential  field  when  enclosed  in  a 
spherical  box  of  large  radius,  a,  on  which  the  wave  function 
must  vanish.   Since  the  differential  equation  is  unchanged, 
all  those  solutions  of  the  scattering  problem  will  satisfy 
the  bound  state  problem  for  which 


(1.3)     sw(|^  a  ^^[e)).o. 

It  is  clear  that  if  we  have  solved  the  scattering 
problem  and  determined  the  function  vi  [Sj  we  can  find  the 
energy  levels,  ^A.>  of   the  bound  particle  by  solving  (1«3) 
for  E;  and  conversely  if  the  spectrum,  H* ,  of  the  bound 
particle  has  been  determined  the  phase  shifts  v^  (^x)  ^^^ 
be  found  from  (1.3). 

These  considerations  lead  us  to  expect  that  the 
Rayleigh-Ritz  method,  which  has  been  so  valuable  for 
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calculations  of  energy  levels  of  bound  systems,  will  apply, 
with  certain  modifications  to  the  calculation  of  phase  shifts 
and  scattering  amplitudes. 


I»2»   Three  Related  Treatments  of  One  Dimensional  Problems. 
The  Collision  Determinant.   Numerical  Illustrations. 

In  this  section  v;e  shall  illustrate  our  methods  by 
treating  the  scattering  of  a  particle  by  a  fixed,  short 
range  scattering  centre.   We  shall  show  that  this  problem 
can  be  attacked  in  three  closely  related  v/ays. 

Let  us  first  simplify  the  Schroedinger  equation  (1.1) 
by  introducing  a  convenient  unit  of  length, -?\(;  ,  and  writing 


(2.1)        OC  ^    — -        , 

K\ 

^^ 

It  then  becomes 

(2.2)     [^        .    ,>_ 

\/(-r 

^1    u(. 

c^«o 

It  will  be  our  problem  to  solve  (1.5)  subject  to  the  conditions 
that 


(2.3 


)     vLio)^Oj        ^tx)  — 3  ^Uvcx  -vt^a^^^^x 


O^      'X.      — ^  <=^     . 

In  particular  v/e    are   interested  in   the   phase   shift  which 
determines    the   scattering  cross   section. 
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Method  1:   The  Energy  as  Proper  Value ♦ 

VYe  know  that  (2»2)  and  (2.3)  have  a  solution  for  every 
positive  VC  •   However,  instead  of  looking  for  the  phase  shift 
corresponding  to  a  given  K.  ,  v/e  can  reverse  the  problem  and 
ask  for  the   K.  corresponding  to  a  given  logarithmic 
derivative,  L,  of  u  at  some  large  distance  a*   This  means  we 
look  for  the  proper  value  \<»    belonging  to  the  differential 
equation  (2.2)  and  the  boundary  conditions 


(2.4)    U.[o)  ^0  )   k  t^)  ^  ^  ^t^). 


) 

This  is  a  proper  value  problem  of  a  standard  type     and  it 
is  easily  verified  that   y<   as  calculated  from  the  equation 

(2.5)  \    l^-u'^+i^^V^Vu?-^  cUc  4  Lu^-^a)  ^o 

is  stationary  with  respect  to  all  variations  of  u  which 
satisfy  the  condition  <Ju[o)cO. 

The  Hayleigh-Ritz  method  is  directly  applicable  to  (2.5). 
y/e  simply  try  for  u  a  linear  combination  of  trial  functions 

(2.6)   U^  ^i  ^lU)  -+  C^U^C:^)  -V C^U^[x) 

(  with  all  U.^  {j>)   -0  )  and  substitute  this  expression  into 
(2.5)»   This  leads  to  an  equation  of  the  form. 


(2.7)  Q   -.X-     [^--\   -^  ^-\  ><"')  ^-  C-  -  0 


where  A;\   ft;',  are  integrals  involving  Ujj^  CUvO.  U.'v 

34«  R»  Courant  and  D.  Hilbert,  Methoden  d«  Math.  Physlk, 
vol.  II,  508,  Berlin,  Springer,  1937. 
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The  stationary  property  of  K   with  respect  to  variations 
of  the  c  »  leads  to  the  following  system  of  equations  for  the 
determination  of  these  coefficients: 


-?  Q 


^  i1(A;j    46.-J   li'-^Cj^o  x.ijt^,,,n. 


(2.8)  ^  Q,  .^ 


As  in  the  case  of  the  bound  state  problem  this  system  is 
compatible  only  if  K  has  certain  proper  values  )i\j^  which 
can  be  found  by  the  solution  of  the  corresponding  secular 
equation.   They  are  approximations  to  the  correct  proper  values 
of  (2.2)  under  the  boundary  conditions  (2.4) •   Their  errors 
are  of  the  order  of  the  square  of  the  errors  of  the  trial 
functions. 

Now  the  phase  shifts  of  the  exact  solutions  corresponding 
to  the  given  L  satisfy  the  equation 

(2.9)      K  cA[\^<l  -^   li^))  ^^    ' 

If  we  substitute  our  approximations  \Lg    in  this  equation  we 
obtain,  to  the  same  (  namely  the  second  )  order  of  accuracy, 
values  for  h  (  i^jt.\  corresponding  to  the  energies   \ij^ 
Then,  either  by  varying  the  pre-assigned  value  of  L  or  the 
"cut-off"  distance  a,  we  can  obtain  n  for  all  other  energies. 
In  this  form  of  the  variational  method,  we  have  followed 
very  closely  the  analogy  betv/een  bound  and  unbound  systems 
pointed  out  by  Pleisenberg.   However  the  quantization  of  the 
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energy  which  entails  the  necessity  for  solving  a  secular 
equation  of  the  usual  kind  seems  artificial  and  can  indeed 
be  avoided* 

Method  2o   The  Logarlthinlc  Dorivatlve  as  Extremum* 

Since  for  a  given  L,   K^  as  calculated  from  (2«5)  is 

stationary,  so  also  must  L  be  stationary  v/hen  calculated  with 

I. 
a  given   vC  •   Using  the  trial  functions  (2.6)  we  are  led 

again  to  (2.7),  but  now  we  take   K    as  fixed  and  L  as 

adjustable  so  as  to  make  the  ensuing  system  of  equations 

compatible*   Therefore  we  exhibit  L  explicltely  and  write 

Q  in  the  form     ^ 

(2.10)  3:., 


-t  -  -  -.   C^  u. 


Vi  "^K 


C^)\'  - 


It  Is  convenient  to  transform  the  coefficients  c  -  to  a  new 
set  d  '   as  follows: 

(2.11) 

In  terms  of  the  d^  ,  Q  has  the  form 

(2.12)    Q^  1    D^- C«.^)    d;  ^i -t  u  a^'' =0 

3--' 

The  stationary  property  of  L  with  respect  to  variations  of  the 
d  \   now  leads,  to  the  equations 

A 


(2a3) 


19 


V\ 


.■^    ^^■^  U'')  ■t'^y-v^-  tA.\  =«5 


k\ 


X.       ^;j    [VC^}    0I3     3  0        ^ 


/I..   -L.^,   _. 


v> 


To  make  the  system  of  equations  (2.13)  compatible  L  must 
satisfy  the  dete.rmlnantal  equation 


(2.14)    A  = 


t>K.      P---- 


I 
t 


D 


Vvlv 


i  0 


Hence  L  Is  a  quotient  of  two  determinants,  each  of  which  Is 

ij  V. 
a  function  of     vs.    • 

Let  us  note  that   A  is  symmetrical  and  of  degree  1 
in  L,  in  contra-distinction  to  the  secular  determinants  in  prob- 
lems of  the  bound  state,  whose,  degree  in  the  unknown  energy 
is  as  high  as  the  number  of  trial  functions.   We  shall  see 
later.  (  sections  I.6.  -  I.7.  )  that  the  structure  of  A-   or 
of  the  corresponding  determinant  occurring  in  method  3  reflects 
the  type  of  the  collision  and  for  this  reason  call  it  a 


collision  determinant 


3r 


•      Its  degree   in  L  is   of  course 


^^*      It  is   related  to  VVheeler'a  Fredholm  determinant,    ref,  ' 
18,    eq.  (45), 
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Invariant  under  the  transformation  (2.11)  which  was  introduced 
only  for  reasons  of  greater  clarity©  Indeed  one  can  easily 
verify  that  the  determinantal  equation  for  L,  as  deduced 
directly  from  (2.10), 

(8.15)   I  C;-^  -V  L  u.;  e^)  ^j  l^)  I  -0 

is  also  of  the  first  degree  in  L,  since  the  matrix 

\l  ^>jC^)  ^'\   t^-)  IV  °^  ^^®  coefficients  of  L  has  rank  !• 

There  remains  as  an  unsatisfactory  feature  of  this 
treatment  the  dependence  of  the  phase  shift,  deduced  from 
^^    and  L,  on  the  cut-off  distance  a,  even  when  a  is 
Y/ell  beyond  the  range  of  interaction.   This  question  will 
be  discussed  in  more  detail  further  on  in  this  section. 
Meanwhile  we  shall  Indicate  another  method  which  does  not 
have  this  drawback* 

Method  5«   The  Tangent  of  the  Phase  Shift  as  Extremum» 
Let  us  return  to  equation  (2.2)   and  (  writing  ^ 
instead  of  K   )  note  that,  in  the  boui:d  state,  ^ 
calculated  from 

(2a6)      I   ^  ^  ""jZc-"^  ^  "^  ^^""^1  ^  ^  ^  ^ 

0 

Is   stationary  relative   to  all  variations   of  u  with    <iu.(0^»Q 
and      Su^oo)    51     0      •     For* 
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and  both   the   integrand   and  Integrated   terms   vanish. 

In   a  collision,   however,    the  wave   function  has   the    form 
(2«3)    at  infinity.      If  we    restrict   our   trial   functions    to 
those  having   the   correct  asymptotic   form,    although 
possibly  v/rong  phase   shift,    then 


(2ol8)  Su     — ^      S^/Ovavi    . 


^.     ^    V<X     ^      O^     ^    ^    ^ 


Hence 


(2ol9)  It      •7au-d«^--l>t 


so  that 


X^    oO 


3t 


^     14    a  JCcUk  v^ 


(2.20)  S  ["i  -^  K  >to^  tj)    ^  0. 

Since   for   the  correct  wave   function  I   :s  0        #    it  follows 
that   the   equation 

(2.21)  X    -V    IC    ;tav^   vl^       =:     iC  jtOc^Vy 

where   vi  ^   is  the  trial  phase  shift,  is  a  stationary 
expression  for  tan  v\      • 

To  apply  the  Rayleigh-  Ritz  method  it  is  essential  to 
write  this  variational  principle  in  a  homogeneous  form* 
Let  us  therefore  consider  wave  functions  of  the  asymptotic 
form 


56«   L»  Hulthen  (  Extrait,  Dixieme  Congres  des  Mathematlclens 
Soandinaves,  Copenhague,  1946  )  has  used  an  equivalent 
equation  to  find  ►!   by  further  restricting  his  wave 
functions  to  satisfy  the  condition  IiO« 
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(2.22)  U      >      ft    5U     Rx     H:     6    ^^  >^  ^ 

instead  of    (2.3)  •      We    then  find  in    the   same   way  as    above 


(2.23)  I  4     AP>l<c.       f\^\L~\kx^ 


1 


as  a  homogeneous  variational  principle  for  tan  u  .   If  we  use 
a  trial  function  of  the  form 


^  ^    l\   SU    ^OC     4  C^  k  ^  4- C^  k  ^ 


(2.24) 

where 

(2.25)    ^t'-L+Cikj-t   .-_     C^K^   __^  B^^^ 


then  since  tan  ^   multiplies  only  A   ,  substitution  of  (2.24) 
into  (2.23)  will  give  a  quadratic  form  of-  the  same  structure 
as  (2.12)  which  in  turn  will  again  lead  to  a  linear  collision 
determinant  for  tan  ^      •  • 

The  drav/back  of  this  formulation  is  that  the  trial 
functions  must  be  of  the  form  (2.24),  (2.25).   For  practical 
purposes  it  may  therefore  be  preferable  to  re-adopt  a  finite 
cut-off  as  follows:   Let  a  be  a  distance  beyond  which  the 
interaction  is  negligible  and  let  us  use  as  trial  function 

(2.26) 

iL  ^  k^"*^  ^  f^  ^^  ^^    -V  B  c^  vex 

where 


(2.87) 


5=    l^'U^)   ^^<^    -  V^'''''^'^)^^-^ 
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A  and  B  are  defined  so  that  u  is  smooth  across  x  ^a.   Clearly, 
since  u  has  the  form  ^     for  x  >  a,  the  Integrand  of  I  Is  zero 
for  X  ^  a  (  see  (2.16)  )  so  that 

(2.28)    X-  ^  uf'M  i--^  -V  4  ^  V(^X)  \  ^  '  dx. 

When  (2.27)  and  (2.28)  are  substituted  in  (2.23),  v/e  have  a 
variational  principle  involving  only   tt*  [x^  ,  whose  asymptotic 
form  no  longer  concerns  us.   The  price  paid,  of  course,  is  the 
reappearance  of  the  cut-off  distance  a. 


The  Well-Depth  as  Proper  Value. -if- 

Before  going  on  to  a  discussion  of  these  methods  let  us 
note  that  if  v^e  write   V(^x)  t:   -T^ -f (x)  >  where  f(x)  determines 
the  shape  and  'X   the  depth  of  the  potential  well,  we  can 
regard  ^  as  the  proper  value  in  the  Schroedinger  equation 

'j\,  »  as  calculated  from  the  equation. 


oO 


cC- 


(2.30)  \      uty)  \  ^^  H  \<''  +  '>-  {i^}    \   UC^)^0 

is   stationary  If 

(2.31)  iu[o)    ^0     ^        Jv.(^)^0. 


"»>  This  idea  was  used  in  connection  with  bound  state  problems 
by  Svartholm  (  ref •  28  ) • 
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Both  K  and  W\  are  now  regarded  as  known  while  the  strength  of 
the  Interaction  Is  the  unknown.   This  is  in  contrast  to 
methods  1-3  in  which  the  interaction  was  regarded  as  given 
and  yL   and  Vi  or  K  and  L  were  determined  as  functions  of  one 
another. 

While  this  approach  has  no  advantages  In  connection 
with  the  Rayleigh-Ritz  method,  it  provides  a  useful  view- 
point in  the  variation  iteration  method  to  be  discussed 
later  (  section  U.S.). 

Connection  between  Methods  1-3; 

Methods  1  and  2   differ  only  with  respect  to  which  of 
the  quantities  L  and  K   is  presumed  known  and  which  is 
unknown.   V/ith  the  same  trial  function  and  the  same  a  they 
must  lead  to  the  same  relation  between  tan  h  and  IC  . 

On  the  other  hand  method  3  differs  slightly  from  1  and 
2.   Let  us  trace,  the  connection  between  the  cut-off  independent 
form  of  method  3  and  method  2: 

In  the  latter,  L  is  determined  from  (2.5),  which  after 
an  integration  by  parts  gives 

(2.32)  -   ^    Zc\  ^       4-  X  H-  L  i-^/O  -  0 

where  I  is  defined  by(2.28).   If  u  has  the  asymptotic  form 
(2. 22) we  can  evaluate  the  two  boundary  terms,  which  leads  to 
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C^  M^  ...  ' 

(2.33)  £.^ 

where   n  f  titCU****  T  B/aN  .   Multiplying  (^.23)  by r+     ^ 

we  obtain 

At   the   particular  cut-off  distances   where  *-  I   ^ 

(2.35)       Ka  ^    '^-^ili-ir, 

(2.34)  becomes 

(2.36)  X     -i-    A^i<.  ytoAA  VI  jj.      -    A    KtctwVl. 

This  agrees  precisely  with  the  result  (2.23)  of  method  3. 

We  now  ask  how  vi  as  determined  by  method  2  varies  with 
a  for  a  given  trial  function  u  which  has  the  correct  asymptotic 
form  but  incorrect  phase  shift.   In  (2.34),  let  us  use  the 
abbreviations 

(2.37)  ^ 

By  some  elementary  manipulations  v/e  can  then  reduce  it  to 

(2.38)  u^   <Svi^  ^^    "      I  ^   • 

For   the  correct  wave   function,    J  -  0,    so  that  by   (2.38) 


f 


«Ai 


Kiev* 
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5  >j  -^  O      ^  except  where  cot  z  =  +  c>o  ,    that  is  where 
\^0\,-\y\j^   "=-  ^"^   •  ^^   those  points  ^y\    Is  undefined.   This 

Is  understandable  from  (2.5)  since  at  those  points  u  r  0  so 

that  L  is  not  determined. 

The  variation  of  dvx    with  z  is  shown  for  different 

values  of  J  in  Graph  1.   For  small  values  of  J,  cot  z  is 

negligible  as  compared  to  l/J,  except  at  the  danger  points 
^^Q^Mx,-  o  =.  yC\\        *    3°  that  in  the  limit  as  J  — >  0,   6n 

becomes  independent  of  a. 

Bilinear  Formulation. 

The  Schroedinger  variational  principle  for  bound  systems 
is  usually  stated  in  either  of  the  tv/o  forms: 


(2.39) 


(2.40) 


where    all   trial   functions    satisfy 

(2.41)      (?u[o'):^0  ^      Su(ob)  :^      0    " 

Let  us  recall  that  a  modification  of  (2.39)  for  the 
case  of  a  collision  led  to  methods  1  and  2,  while  (2.40) 
was  the  starting  point  for  method  3. 
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We   can  replace  (2, 39)  and  (2.40)  by  the  following 
equivalent  statements,  involving  bilinear  instead  of 
quadratic  forms : 


or 


ji^ 


'    \\^     €^   :s  0 


(2.43)  ^^         ^^ 

where  all  trial  functions  U,  and  U^  satisfy 

v      ^ 

(2.44)  <?u4t))=0^     5ui.(^)   :.  O  ;  =  \^i. 

0  in  (2.42)  and  (2.43)  indicates  independent  variations  of 
\X^    and  U^.   The  truth  of  (2.42)  and  (2.43)  is  directly 
verified  from  the  Schroedinger  equation  and  the  boundary 
conditions  (2.44) • 

Clearly  the  corresponding  variational  principles  for 
collision  problems  can  be  similarly  modified.   Without 
going  into  details  we  shall  merely  state  the  results. 

Methods  1  and  2: 

if  the  trial  functions  satisfy 


(2,46)     5l^;(^)30 
(  Compare  to  (2,5).) 
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,    -  -  ^;  ^ 


Method  3: 


(2.47) 


If  the  trial  functions  satisfy 

(2.49)  '  }    A^     \;,2.. 

vt^ ^    A;  s>'—  »<- X   -I  6^-  (i^  v<  dfJ 

(  Compare  to  (2.23),  (2.24).)   It  Is  not  at  once  apparent 
that  tan  Vj   as  calculated  from  (2.47)  is  symmetrical  In   U.  ^ 
and  Iv  .   That  this  is  so,  however,  follows  from  the  fact  that 


'v-idO 


-  ^ai-  f\^\s^.  h^%^\L  _ 


In  one  dimensional  problems  there  is  no  advantage  in 
using  the  Rayleigh-Ritz  method  in  conjunction  with  the 
bilinear  variational  principles  instead  of  the  quadratic 
ones.   But  wo  shall  find  the  bilinear  version  useful  in  three 
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dimensional  problems  (  sections  1.4,  -  1.5.  ),  where  the 
Schroedinger  equation  has  more  than  one  solution,  and  in 
connection  with  the  variation  iteration  method  (  sections 
II. 8.  -  11.10.  ).  . 

Numerical  Illustrations ♦ 

Using  the  m.ethods  2  and  3  we  have  calculated  ^  S 
scattering  of  neutrons  by  protons  at  zero  energy.   We  have 
assumed  the  interaction 

(2.51) 

(  X  measured  in  units  of  A^^Z'30  y^io'        cm  ) ,  v/hich  gives 
the  correct  binding  energy  for  the  deuteron.   The  wave 
equation,  in  this  case,  is  simply 

(2.52)  i^      ^    ^(:,)  i^  ^0. 

and  can,  of  course,  be  solved  exactly.   For  x  >  1  the 
solution  has  the  form  u  s  A  (  x  +  X  )  where  X  is  the  desired 
limiting  ratio  of  tan  yj  /  ^    ,    as  K  approaches  zero. 


As  trial  functions  v/e  have  used 


(2.53)  a»,  C  l  ^g^ 


I 

with  n  s- 1,  2  and  3.   Since  (2.53)  does  not  have  the  correct 
asymptotic  form,  it  was  necessary  to  use  method  3  with  a  cut- 
off, a*   The  interaction  vanishes  at  x  s=.  1,  which  is  the  most 
appropriate  choice  for  a.   However  to  study  the  dependence  of 
X  on  the  cut-off  distance,  we  have  made  calculations  also  for 
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larger  values  of  a.   Clearly  a  function  of  the  form  (2,53) 
can  approximate  the  correct  wave  function  better  in  a  small 
than  in  a  large  interval,  so  that  as  a  increases  we  expect 
our  approximation  for  X  to  get  worse. 

The  results  are  listed  in  Table  I: 

Table  I, 
Approximations  to  X  g  -^^  "^^^^  ^  /  K  . 


C/JcO/cJt 


s;        -    2.080 


a 

^1 

X 

C  X  +  C  X 

C  X  +  c 

1              ^^ 

X-       3 

X   +  C  X 

1 

Method  2 

Method  3 

Method  2 

Method  3 

Method  2 

Method  3 

1 

-  3.957 

+  1.338 

-  2.112 

-  2.126 

-  2,083 

-  2.084 

2 

+  4.043 

+  1.338 

-  2.853 

-  4.039 

-  2.201 

-  2.206 

3 

+  2.416 

+  1.338 

-  6.593 

-  51.606 

-  2.385 

-  2.520 

4 

+  2.011 

+  1.338 

+968.870 

+  7.634 

-  2.936 

-  3.315 

These  results  exhibit  a  number  of  significant  features: 

1.  The  methods  2  and  3  lead  to  errors  of  the  same 
order  of  magnitude.   (  In  the  present  example  method  2  is 
somewhat  better*  ) 

2.  As  expected,  the  results  grow  progressively  worse 
with  increasing  a. 

3.  The  dependence  on  a  becomes  weaker  the  larger  the 


31 


number  of  parameters  (  except  for  the  1.338  coluran;  see  5. 
below  ) • 

4.  There  appear  singularities,  owing  to  the  wrong 
asymptotic  behaviour  of  the  trial  functions.   As  the  number 
of  parameters  increases  the  singularities  move  outwards,  to 
infinity. 

5.  The  trial  function  ex,  in  method  3,  leads  to  a 
result  independent  of  a  but  of  the  wrong  sign.   This  can  be 
understood  as  follows :  c.  x  is  the  correct  solution  corresponding 
to  no  interaction.   It  has  (  accidentally  )  the  correct 
asymptotic  form  (  though  wrong  X  )  and  therefore,  by  the  theory 
of  method  3,  must  lead  to  results  independent  of  a.   In  fact 

v/e  have       ^ 

(2.54)        "        , 

all. of  which  are  independent  of  a.   Using  these  in  (2.23)  we 
obtain  for  the  determination  of  X  c  tan  n  /  \^ , 


(2.55)    -    )    V(^)  Ui^'tx)  "^   -     c,^  J::^^  i|/ 


which  is  precisely  the  Born  approximation  for  tan  M  .   For 
small  V  the  latter  gives  good  results.   But  as  V  increases 
(  tan  V)  )g()jt4.   has  a  fixed  sign  while  the  true  tan  >n 
eventually  changes  sign.   This  has  occurred  in  our  case, 
hence  the  v;rong  sign  in  the  1.338  -  column. 

In  summary  we  may  say: 

(a)   Method  1  is  the  least  natural  to  collision  problems. 
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(b)  Method  2  appears  to  give  slightly  "better  results 
than  method  3;  but  it  suffers  the  aesthetic  flaw  of  necessarily 
involving  a  cut-off  distance. 

(c)  Method  3  is  closely  related  to  method  2.   It  does 
not  necessarily  involve  a  cut-off  and  yields  the  Born 
approximation  for  the  phase  shift  as  a  special  case. 

(d)  The  Rayleigh-Ritz  principle  applied   to  either 
method  2  or  3  leads  to' a  collision  determinant  instead  of  the 
usual  secular  determinant. 


1.3.   Variation  of  the  Phase  Shift  with  Energy.   A  Stationary 
Expression  for  the  Width  of  a  Resonance  Level. 

We  have  seen  that  methods  2  and  3  lead  to  stationary 
expressions  for  the  phase  shifts  corresponding  to  arbitrary 
energies.   It  follows  that  it  must  also  be  possible  to 
obtain  a  stationary  expression  for  ^^  I  JL  v^         which  determines 
the  variation  of  the  scattering  cross  section  with  energy. 

It  is  indifferent  which  of  the  tv/o  methods  v/e  use,  but 
method  2  seems  more  convenient.   We  introduce  the  definition 

(3.1)  n^[u)     E  ^^  [-    ^T   ^    ■^''-''  -A^')-]^   ^W^^IO), 

where    Lu      is    the    logarithmic   derivative   of   the   correct 
function    U, J.   belonging   to   the   energy    K    •      By  the   Schroedinger 
equation    (2.2),  H  »   f  U.  jv ^  zO  $    ^^^  ^7   the   stationary  property. 
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(2#5),   we   have   rigourously 

(3.2)    M^^(u)  ^    MkC^^-k)   -   oCSu";)^ 

for  any  function  u.  -  w.^.  -t  Su.^^  •   Differentiating  (3.2)  with 
respect  to  K  we  obtain,  to  the  first  order,  in  S^y^  $ 

(3.3)  __JllJ^   ti  0 

Hence   by    (3.1),    we   find  the    equation 

(3.4)  2k    \      ^"^         ^       ^     ^"^^)     "^ 

o 

which  is  a   stationary  expression  for  di'\^   /  t?(  K 

From  this  equation  we  can  determine  ^^J e{(\^)a3    follows. 
We  note  that 


(3 


.5)  L,^  ^  ^cA  [yc^^Y^ 


Hence 


Al^    ^    ^  [\i^-v  M  k)  -- c-^'^^'^(^v4<i  +  y)[^Ka -+K.  Jl!5\ 


(3.7) 


so  that 

v./  V        V       V  y     — — — 

a  ^ 
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We  may  observe  that  this  stationary  expression  for 
involves  only  the  wave  function  and  the  correct  phase  shift 
V] ,/  ,   while  all  reference  to  the  interaction  has  disappeared. 

Let  U3  check  that  if  <Av*j^^dH»^Jls  computed  from  {3»7), 
using  the  correct  v/ave  function  i^y.  9    ^^^  apparent  dependence 
of  c^l^M*^ jon  a  disappears,  as  it  must,  as  long  as  a  exceeds 
the  range  of  interaction,  ^^. 

For  X  >  Q.^  ,  the  correct  H^  has  the  form 

(3.8)  Uj^[x)  ^  C  SU^KaH^vt.) 


Hence 


a 


Differentiating  with  respect  to  a  we  find 


^v<.       ^K. 


(3.10) 


^  I  c^  il[K/d-v^K)]  ^;s»/-^(Kaiy)  -ij 


which  proves  that  <li</cl(^v\  is  independent  of  a. 

As  an  application  of  (3»7)  we  shall  derive  a  stationary 
expression  for  the  level  width  at  resonance. 

In  the  neighbourhood  of  a  resonance  energy,  defined  by 
the  condition  that  »1  (^  ^  lv\  4  — ^iT  >  tan  vi  has  a  singularity 
of  the  form 
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(3.11) 


UJUa 


where  K;|  is  the  resonance  energy.  \         l3  called  the  level 
width  because  the  cross  section  in  the  vicinity  of  K^  is 
given  by 

(3.12)   .        0-  ^      Jb 

Now  near     K^ 

(3-13)       nx  M^-t-^V    +  -J^    ( '^''- '^^''^j  • 


Therefore 


(3.14)  -tOw  M  K    ^     — 


d-^v. 


K[0     . 


1-  ./     V 


k^-'K 


Hence,  comparing  (3.11)  and  (3«14)  and  using  (3.9)  with  v\   ts,[vv-t— \  il 
we  find  the  following  stationary  expression  for  •!•  i   : 


ir  =    - 


(3.15) 


V 


<^i^i) 


U-K 


4k*         ^ 
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I«4«   A  Variational  Principle  for  the  Scattering  Amplitude. 

When  dealing  with  a  wave  function  of  given  angular 
momentum,  a  correspondence  to  a  bound  state  problem  can  be 
established  because  we  know,  a  priori,  that  asymptotically 
the  solution  has  a  radial  logarithmic  derivative  independent 
of  the  angle.   But  if  we  want  to  treat  the  whole  incident 
plane  v/ave  at  once  no  such  information  is  available.   We 
have   already  seen  that  only  method  3  need  not  involve  an 
enclosure  and  that  the  Born  approximation  for  individual 
phases  can  be  deduced  from  it.   Both  facts  suggest  that  it 
v/ould  be  the  most  promising  attack  on  the  scattering  amplitudes 
of  three  dimensional  problems. 

The  wave  equation  for  the  relative  motion  of  two  colliding 
particles  of  equal  mass  M  is 

(4.1)  1-1\%  [)(\<)]  ^[\<)  .  e  ^(^), 

where  R  is  their  relative  position  vector  and  U  (R)  is  the 
interaction  potential.   Setting 


■^o 


2^ 


t> 


(4*2) 


vvW-*(^.'), 


where  ^,  Is  our  unit  of  length,  it  becomes 
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(4.3)  U'^^  \<^   -  V/(irJ^  v^(r)  =^0 


We  shall  nov/  follow  a  procedure  corresponding  to  the  bilinear 
form  of  method  3  (  see  (2.47)  -  (2.49)  ).   We  define  O^^^    ,    and 
vJU  ^^    as  the  two  solutions  of  (4.3)  which  correspond  to  plane 
waves  incident  along  -K  and  \<^  and  have  the  asymptotic  forms 

(4.4) 

Here  ^Kj  ^  -  [U^j  c[i^)  and    ^[*^-:.^)  ^^  ^^®  scattering  amplitude 
in  the  direction  1<  of  a  wave  incident  along  K^  •   Corresponding 
to  (2.47)  we  now  define 

In  virtue  of  (4.3),  we  find  after  an  integration  by  parts, 

s 

where  S  is  a  large  sphere,  and  ---  denotes  differentiation 
along  the  outward  normal. 

Now  let  us  admit  only  trial  functions  of  the  asymptotic 
form  (4.4),  but  v/ith  possibly  incorreat  f  *3,  so  that 
asymptotically 


(4.7)    ^^K^   — ^   ^^t'^^C^^.,  V<) 


^ 
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v/here     o  \      \^Li^)         ^^    ^^®  error  In  the  scattering 
amplitude  of  the  trial  function  vV^  •   Clearly  only  the  plane 
wave  part  of  \J^^    contributes  to  the  surface  Integral  (4«6) 
3  o   that 


SI'..,  -\r'''k   '- 


e. 


"^1  <:  9^''>f     \ 

Dlrac    has  shov/n  that  the  factor  multiplying  d  A     (Kj^KJhas 
a   8  -  function  like  character  for  large  values  of  /t  ,  so 
that  the  Integration  gives  simply 


(4.9)  SIk,,k,  --Mir  M''^(Km'^')- 

Since  for  the  correct   Mj'k^  *   -^  >^ii  t<   ~  ^   *    ^^   follows 
that  ^{^2.     U*\\  9    as  calculated  from  the  equation 

(4.10)  X^^^^,   +  Kir  \^^  C>S'^i)    -  4Tr-{Cs*^0 


is  stationary  with  respect  to  independent  variations  of  ^'-Ki 
and   d/^  ,  subject  to  the  conditions  (4«7), 

Symmetry  Properties* 

Equation  (4»10)  involves  two  trial  function,    ^'-k  ^^^ 
vJ/y(^   .   These  differ  by  a  rotation  in  the  coordinate  space 

37.   P.  A.  M.  Dlrac,  The  Principles  of  Quantum  Mechanics, 
(  Clarendon  Press,  Oxford  1947  ),  third  od,,  p.  191. 
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and,  in  general,  also  in  their  functional  form  (  to  the 
first  order  )•   Let  us  indicate  the  functional  form  by  super- 
scripts  and  write  explicitly   H-w*  ^  ^t,   ^*^^  ^^^  trial 
functions,  and  J^    (l<,,^Kv)f'or  the  f  obtained  from  (4.10)  by 
using  vV  j^  first  and  O/^^    second  in  (4«5)«   Integrations 
by  parts,  analogous  to  (2.50),  reveal  t}ie  following  symmetry 
properties  which  are  not  at  once  apparent  from  (4.10): 

(4.11)   t^^i^.,K,]     --  f''\^^.,^.)    , 


(4.12)    f-U-'^^,'*0  -  {"^C-^MV^^j. 

(4.11)  implies  that  f  is  invariant  under  interchange  of  the 
functional  forms  of  4^-^.  »    ^nd   *^Ki  »  ^^  will  be  verified 
below  in  the  case  of  the  Born  approximations.   (4.12) 
guarantees  that  the  scattering  kernel  -^(-Kj..\<^^  is  symmetrical 
when  computed  by  (4. 10). 

Connection  with  Born  Approximations. 

The  simplest  admissible  trial  functions  are 


(4.13)  ' 

Mr  Kj.  ^  e 


38.   Ref.  18,  eq.  (2.52) 
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Substituting  these  in  (4.10)  we  obtain 

(4.14)  -^e"^^'"'  V  e'''^-'    .    Mn  {i^,^.) 

which  is  identical  with  the  result  of  the  first  Born 
approximation* 

V/e  obtain  the  second  Born  approximation  if  we"  use  one 
first  Born  approximation  trial  function  and  one  plane  wave; 
i«e»  either 


or 

(4.16) 

in  agreement  with  (4«11)«   Higher  Born  approximations  are 
obtained  in  a  similar  way. 

Thus,  apart  from  re-deducing  the  Born  approximations 
from  a  variational  principle,  nothing  new  is  gained  by  using  theso 
iterated  functions  in  (4.10) . 

Application  of  the  Raylei/?,h-Ritz  Method* 

We  can  use  the  Rayleigh-Ritz  Method  in  conjunction 
with  (4.10)  to  determine  the  scattering  amplitude  >(  [k^  K^V 
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This  Is  analogous  to  method  3  (  section  1.2.)  for  calculating 
the  tangent  of  the  phase  shift.   For  simplicity,  we  use  the 
same  functional  form  for  ^^u    ^nd  ^^     and  write 

(4.17)  ^    *^l  ,^  "^         ^\  i^  / 

where  the  trial,  functions  ^\c\^      have  the  asymptotic  form 


(4.18)  Ul^r) ^  i;  U    k') 


Clearly  for  use  with  trial  functions  of  the  form  (4.17)  instead 
of  (4.4),  (4.10)  must  be  replaced  by  the  equivalent  homogeneous 
equation 

(4.19)  1r^^^, -f  M-M  c,^   ^^{^[^^,^.)   -  i|ir  C  ,''  f  (  K,  , /i.)  ^ 

where  J-H..  vt,  Is  defined  by  (4.5). 

When  (4.19)  is  written  out  in  terms  of  the  trial  functions 
it  assumes  the  form 

(4.20)  "l.     Ar.  c-c-    .  4u   {(k,,k,)  c,^ 

» 

where  we  may  take  "^^  *i  ^  ")/»  •   The  stationary  property  of  f 
relative  to  variations  of  the  C*  then  leads  to  the  equations 

Z.   A  I)  c-    -  M«   J»/k,    ti)  c. 


^  A;j  ct  -  p        £  «4j..- 
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They  are  compatible  only  if  the  determinant  of  the  coefficients 
vanishes,  that  is  if 

f»n-Hr;li{^,«.)  A,v  -----  Am 

^U  ftiv  Alt, 


(4.22) 


^{^t,^)   £ 


"Kl 


I 
Ah^ 6 


KU 


-  0 


f(  ^i,\L.   )    is  uniquely  defined  by  (4.22).   The  fact  that  it  is, 
in  general,  complex  is  reflected  by  the  collision  determinant 
A  f  K  \L  \     whose  matrix  is  symmetrical  but  not  hermitian. 

On  varying  the  angle  between  -^y    and  k^,  one  can  obtain 
stationary  expressions  for  the  scattering  amplitudes,  and 
hence  for  the  differential  cross  section,  in  all  directions. 


I«5«   Non-Central  Forces* 


It  is  well  known  that  if  two  particles  interact  by  means 
of  a  potential  which  depends  only  on  their  distance  of 
separation,  the  relative  angular  momentum  is  conserved.   This 
fact  is  the  basis  of  the  phase  method  in  collision  problems, 
v/hose  variational  formulation  has  been  discussed  in  sections 
1.2.  and  1.3. 

However,  the  hypothesis  that  nuclear  forces  are  purely 
central  is  untenable  in  view  of  the  quadrupole  moment  of  the 
deuteron,  and  has  been  replaced  by  the  assumption  that  they 
are  a  mixture  of  central  and  non-central  forces.  Rarlta  and 
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Schwinger    have  worked  out  the  combinatorial  aspects 
introduced  into  the  2-body  problem  by  the  existence  of  such 
non-central,  so-called  tensor  forces.   In  particular  they 
have  shown  that  the  angular  momentum  is  no  longer  conserved, 
and  that  the  eigenstates  of  the  Hamiltonlan,  instead  of  being 
S,  P,  D, . . . .  states  are  now,  in  general,  mixtures  of  such 
states;   Most  important,  in  the  neutron-proton  system  at  low 
energies,  is  the  S  +  D  mixture • 

In  our  method  of  section  1.4. ,  the  conservation  of 
angular  momentum  did  not  have  to  be  assumed  and  thus  it 
applies  without  change  to  collisions  involving  tensor  forces. 
This  method  deals  with  all  angular  momenta  at  once.   It  may 
hov/ever  be  desirable  to  treat  each  elgenstate  separately. 
This  can  be  done  in  a  way  related  both  to  the  three  dimensional 
procedure  of  section  1.4.  and  to  either  method  2  or  3  of 
section  1.2. 

As  mentioned  abovd,  it  has  been  shown  from  the  trans- 
formation properties  of  the  Hamiltonian  that  there  are  solutions 
of  the  Schroedinger  equation  consisting  of  mixtures  of  S  and  D 
states.   These  solutions  are  proper  functions  of  the  Hamiltonian, 
H;  the  total  angular  momentum  J;  its  z-component  *^x  ;  the  total 
spin  S;  and  the  parity  operator,  P.   In  the  scattering  of 
neutrons  by  protons  the  incident  plane  wave  includes  both  S 
and  D  waves  in  which  these  operators  have  the  same  proper 
values.   The  amplitudes,  Rs  *  P>>  of  the  outgoing  S  and  D 
waves  are  related  to  the  amplitudes  ^$  #  ^9  ^^  ^^®  Incoming 
waves  by  a  unitary,  symmetrical  scattering  matrix  S.   Thus 
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(5.1)    ft  ^  S'   o( 


where 


'^n^p 


(J 


^ 


(5.2) 


We  know  o^^  and  ^p.   To  find  the  differential  cross  section 
we  must  determine  the  elements   S^*'  v/hich  give  A^  and  P  D 
by  (5.1).   Actually  it  is  sufficient  to  determine  two  of  the 
(  complex  )  elements  S^'v  since  due  to  its  symmetry  and  unitary 
character,  the  matrix  S  involves  only  three  independent 
constants. 

Consider  the  expression 

ON  (V) 

where     vl^    and     v^    are    the   solutions   of   the  Schroedinger 

equation. 


(5.4)     (v%K--V)^-_0^ 


which  have  the  asymptotic  form 


39.  See  ref.  (18). 
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t    — ^ 


D 


(5.5)    .U)         ,   <•»      -(^        .iw  „-'"  A   fN  -^ 


XL, 


Here  -^s  >  *^^9  describe  the  normalized  angular  and  spin 
functions  of  the  S  and  D  v/ave  respectively.   Now  let  us  use, 
in  (5.3),  trial  functions  whose  asymptotic  variation  is 


(5.6) 


Then 


(6.7) 


Hence,  to  the  first  order,  we  have 


(5.8) 


or,   in  the  matrix  notation  of   (5.2) 

(5.9)  ^  "  ^ 


ivl 


s    ^Ji   ,t'"-S-ci. 
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This   enables  us    to   obtain   all  matrix  elements    ^X\    •      For 


s. 


example,  to  find  -^^^   we  use  trial  functions  with  O^p  =>0  ,  . 
S 


ol^  'O     in  which  case  the  right  hand  side  of  (5«9)  becomes 


(5.10,    If      [OC,   0)(^^^   J(,«j  =  ^C.,  0^0  ^.^^ 

and  an  application  of  the  Rayleigh-Rltz  method  will  yield  the 
required  matrix  element. 

Let  us  note  here  that  in  this  way  we  obtain  a  stationary 
result  for  the  differential  a3  well  as  for  the  total  cross 
section  (  The  stationary  character  of  the  differential  cross 
section  is  lost  in  the  method  of  proper  phases,  described  below.) 

The  Structure  of  the  Collision  Determinant* 

Suppose  v/4  want  to  calculate  the  element  ^u  by  the 
Rayleigh-Ritz  method.   We  use  as  trial  functions 


1 


40 •      For  calculating    S^     we  need  not  use   different   functional 
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where  the  radial  functions  U.   and  kr    vanish  rapidly  as 


/t  ^  C3-0  •   Thus  each  trial  function  consists  of  an  S  and 

a  D  part.   Using  the  functions  (5.11)  in  (5.9)  we  find  that 
the  tensor  character  of  the  potential  introduces  elements 
connecting  all  coefficients  of   vK»,  with  those  of   v^^ 

The  unknown   S((   is  introduced  by  the  right  hand  side  of  (5.9) 

(il  {^l 

and  occurs  only  as  factor  of  ©/j  c^<  • 

On  using  the  stationary  property  of   '^U  ,  we  are  led 
to  a  system  of  equations  precisely  as  in  section  1.2.   These 
are  compatible  only  if  the  determinant  of  the  coefficients, 
b^y^         ,  vanishes.   In  the  present  case   Ajy  has  clearly 
the  following,  symmetrical  form: 


0) 

•  %- 

jl. 

-^ 

s 

u.) 

Ik 

s 

_p_ 

/ 

A«^ 


(5.12) 


C^u  = 


.M 


5f-»Vu 


I' 


n 


K 


X 


X 


^•*1^u  X 


'i 


X 


X 


n 


"  I 


\> 


s 


D 


I 


o; 


% 


V 


*0       ,  i^) 

forma  for  ^fc     and  MT^  •  v;e  use  them  here  because  for 
calculating  off-diagonal  elements  of  S,  different  functional 
forms  obviously  are  required. 
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Because  of  its  symmetry  L^x\    -^(Au)   #  where  J^^^     is  the 
large  upper  right  hand  block  of   A^  •  ^u     is  therefore 
determined  by  the  linear  equation  A-^  r^O   •  /S\|   comprises 
two  blocks  of  elements  connecting  S  and  D  functions.   These 
describe  the  interaction  between  the  S  and  D  state.   In  the 
limiting  case  of  vanishing  tensor  forces  they  become  zero 
and  ^        assumes  the  degenerate  form 


(5.13)    Au 


Hence     >,      is    determined  from  the   equation 

S_ 

X 


(5.14) 


;t*3^„ 


S 


'   0. 


whose  elements  no  longer  involve  any  D  functions.   This 
expresses  the  fact  that  the  S  wave  is  now  decoupled  from  the 
D  wave.   In  fact  (5.14)  has  the  same  structure  as  (2.14), 
which  v/as  arrived  at  by  considering  only  the  S  waves. 
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Proper  Values  of  the  Scattering  Matrix. 

(r)      /i; 

The  unitary  matrix  S  has  two  eigenstates   CP   and  (J 

with  the  property  that 


(5.15)    (?)5   S-o^''   -   X""  Oi'*    [^     ^"-"^ «-^y 

The    /v      are    the  proper  values    of   S   and,    "because   of   its   unitary 
character,    may  be   written  as 


(5.16)   X>'^  -e^-'l^  ,    )C^^-  - 


Thus  the  eigenstates  have  the  asymptotic  form 


<^-^')         ^dfi   ^-^    .e""l^  e^--^^ 


D 


'^  -^  1^ 


A 


II 


.  ^-(k;..^^)[c^^,^>J1^^  4^11^3 


(5.18)   c^^""*  ^>.^;jn^  .  sU[^^,^ll)[df'Sis•^Jf^,) 

Cp   and  Cp  are  thus  standing  waves  with  real  radial  functions 
and  for  that  reason  they  are  simpler  to  deal  with  than  other 
states  which  are,  in  general,  complex  linear  combinations  of 
Cp    and  iP^  *  .   Furthermore  Schwinger    has  shown  that  the 

41.   J,  Schwinger,  unpublished  lectures. 
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total  cross  section  of  the  S  +  D  mixture  of  a  plane  wave  can 
be  expressed  in  terms  of  the  two  constants  V)^   and   viy  only, 
namely 

(5.19)    5 \    \^^^      vjx  +  5^!^  yy^    J 

It  is  therefore  of  interest  that  tanh^.  and  tan  ^^  can  be 
directly  obtained  by  our  method.   In  fact  this  problem  is  very 
much  akin  to  that  of  finding  tan  vi     in  the  case  of  central 
forces  {  section  1. 2,  ). 

We  could  proceed  from  the  most  general  equation  (5.9) 
but  the  following  approach  is  slightly  simpler: 

The  correct  CP     and  Cy       both  have  the  form 


(5.20) 


C/>  and .  CP 

^\^,  -±^   ^^,  <r!^^i,u,,,...}il 

where  the  a  .  ,  b  *  ,  u*  and  w*  may  be  taken  as  real  and  ^  -  ti-tc^/)        . 
Consider  now  the  expression 

(5.21)    X   *  W"*"  {_^\  v.""-  V)  v^  (ct  v) 

I 
1 


where  we  admit  only  such  trial -CP5  >  whose  variations  have 


the  asymptotic  form 


(6.22) 
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Proceeding  in  the  usual  way,  one  finds 


-Ul.^l^s 


(5.23) 


"h'i'  ■"if.t 


Therefore,  corresponding  to  (2.23),  we  find 

(5.24)  r+  K[a,^L -«- A-^i)  -^  K (5 v^  +  ^^^J  ;t<u.  h 

as  a  stationary  expression  for  tan  ^Ta;iTr» 

When  we  substitute  a  trial  function  of  the  form  (5.20) 
in  (5.24)  we  obtain  an  equation  of  the  type 


(5.25)   2. 


I    it- 


^^  ^]  -V  Q^j  ^;  ^')^'^;']^;l-\'  ^[a^-tK)^^. 


Proceeding  in  the  usual  way  to  a  set  of  linear  equations  and 
their  compatibility  condition  v/e  find  the  following 


determinantal  equation  for  tan 


(5.26) 


a  - 


r 


Ti.-xiu,,,    ?,^ IQ,, 


i\. 


p 


I  Qi. 


o» 


Q 


tl 'X^~Kt<U.Wj   tSti  .  , 


Qx. 


I  ' 

!  \ 


% 
I 

f 


s 


^  0 
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Clearly  this  is  a  quadratic  equation  for  tan  n  ,  yielding 
approximations  for  tann^  and  tan  n—  •   The  elements  Q;^\ 
arise  from  the  interaction  of  S  and  D  functions.   In  the  absence 
of  tensor  forces  they  vanish  so  that  A   factors  and  we  are 
led  back  to  tv/o  independent  problems  for  the  S  and  D  wave 
(  compare  (2#14)  )• 

From  the  two  solutions  for  tan  n   we  can  find  the  total 
cross  section  by  (5.18),  whose  error  will  be  of  the  second 
order  if  our  trial  wave  fxinction  has  a  first  order  error.   We 
can  also  find  the  differential  cross  section  to  the  same  order 
of  accuracy  as  our  trial  function  by  solving  the  linear  equations 
derived  from  (5.25)  for  the  d;   and   ^a  ,  thus  obtaining  the 
wave  functions  CP       and   C^  '  and  then  taking  that  linear 
combination  whose  incident  wave  agrees  with  the  incident  S  and 
D  waves  contained  in  a  plane  wave.   To  obtain  the  differential 
cross  section  to  the  second  order  of  accuracy,  we  must  fall 
back  on  the  more  general  but  more  complicated  equation  (5»9)» 

r 

Numerical  Illustration. 

By  using  the  method  of  proper  values  of  the  scattering 
matrix,  we  have  calculated  tan  hj-  and  tan  hi"  for  neutron- 
proton  ^Si-f  ^  ^l     scattering  at  2.82  I.ffiV.   We  have  used 
the  constants  of  Rarita  and  Schwinger    to  have-  a  comparison 
with  Schwinger »s  numerical  integration    for  the  same  energy. 
Thus  the  following  interaction  potential  was  used: 

(5.27)     V  =   1^{  \[0  +2."'^'j(r)  Sa^    ^ 
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where 


fH-   -\  ;   ^^ 


I 


(5,28)  •:-  0  A  >i 

and   S|2^  is  the  spin  dependent  tensor  operator. 


V 


(5*29) 

A 


-»3 

r   is  measured  in  units  of  ^q  ^     2.80  x  10    cm. 

Schwinger   has  shown  that  the  Schroedinger  equation 
leads,  in  the  present  case,  to  two  coupled  equations  for  the 
radial  functions  U.  (^)  and  ^[^)   of  the  S  and  D  wave,  viz. 


(5 


•^^^iX-'^'V^^i^'T)  "' 


(5.31) 


where  Jl    and  Q  are  defined  in  (5.28)  and  \\c-   \- 1  *^   .   It 
Is  simpler  to  base  our  variational  procedure  on  these 
equations,  rather  than  to  go  back  to  the  original  three 
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dimensional  Schroedlnger  equation.   We  define 


(5.32) 


'"l(C''''-i)"^>(r-''^" 


H 


If  we  use  trial  functions  with  the  asymptotic  form 


k 

(5.33) 

U; 

we   find 

(5.34) 


Now  the   exact  solutions  which  we   are   looking  for  have    tan    h  ^ 
r    tan  ^  Fs    "^  tan   h  .      Hence    the  equation 


-K[As''-tAp^J  ;^t<u^  vj 


wherd  "t"  Indicates  the  trial  function  as  usual,  gives  a 
stationary  expression  for  tan  vj  • 

We  have  used  as  trial  functions 


55 


^'^--^'^^  K-t,^.3 


(5.  36)         ^       ^  ^  \    ^^\ 


A  >  -^h  >  ^^^^4'  ^^^  ^^^  ^A"  ^''®^®  determined  from  the 
requirement  that  at  A  :=•  1  the  functions  (5«56)  pass  smoothly 
into  the  free  particle  functions 

(4) 


(5.37) 


As  checks  on  our  calculation  we  have  at  each  stage  computed, 
by  the  variational  method,  the  S  and  D  phase  shifts  when  tensor 
forces  are  absent  and  compared  the  results  with  the  exact 
values  which  are  easily  obtained  for  the  uncoupled  states* 
The  results  of  our  calculations  are  collected  in  Table  II. 
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Table  II. 


S  •♦•  D  Scattering  with  Tensor  Forces » 


Uncoupled 

Coupled 

Trial  Functions 

XxXM,h 

tou^  VI  s  10^ 

A^X/vvVjx 

:;fcaM,yji[ 

us     C.  A    +    Cjj^A  "1 

6.9850 

-   0.3199 

6.986 

-    0*0003 

w    z  djA^               ) 

U   i  CjA   +    c,,A      ^ 
w   *  d^A^    +  d^  A**  J 

6.9850 

-   0.3070 

-   1.525 

-   0.0002 

1.              I 

U    ■&    Cj  A    +    C-j^  A     +    CjA-j 

w  c   d  A,     +  d^A              J 

• 
9.1398 

-   0.3070 

-    1.425 

-   0.0002 

Correct  Functions 

9«2314 

-   0.3069 

-    1.435 

0.0000 

It  is  seen  that  the  last  pair  of  trial  functions  gives 
good   agreement  with  Schwinger's  numerical  integration. 

1.6.   Neutron  Deuteron  Scattering. 


So  far  we  have  been  concerned  only  with  two  particle 
scattering  processes..  We  shall  show  that  the  methods  which 
we  have  developed  can  be  extended  to  many  particle  collisions. 
In  this  section  the  simplest  collision  of  this  type,  namely  that 
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between  a  neutron  and  deuteron,  will  be  considered*   Our 
procedure  will  be  modeled  after  method  3  in  section  I. 2., 
but  method  2  is  equally  well  suited  for  many  particle 
collisions. 

The  wave  equation  for  the  system  of  the  two  neutrons, 
"1"  and  "2",  and  the  proton,  "3",  is  by  (0.2)  and  (0.5) 

(6.1)  l^T  +  \fl^^)^M[^J)   ^  Visi)}  ^   ^^v.-t^^)^ 

where  T  can  be  written  as 


(0.5)  T-^^  JL  f  ^V-^4  ^//) 


HI 


To  remove  the  cumbersome  factor  in  T,  we  make  the  follov/ing 
definitions: 


IV 


(6.2) 


% 


^-D/ 


4ii 


r 

Eh 


42«  For  the  definition  of  the  relative  coordinates,  see  (0.3) 
and  (0.4). 
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Thus 

(6.3)       D.  -  (Vr;4^  V) 


We  now  define 


(6.4)  ?  E   -  D  -  ^.vU^^-Uj,  ^si-^^^ 
SO  that  16.1)  can  be  written  In  the  form 

(6.5)  I'^^O 

Now  consider  the  expression 

(6.6,    X,  ^^    ^'?^  [lh)[.^^.], 

where  V^  Is  the  volume  bounded  by  the  five  dimensional 
surfaces. 


(6.7) 

y  (V) 


Clearly 
(6.8) 


We  aim  to  convert  this  Into  a  surface  Integral  over  the  surfaces 
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2.      and  2.a   ^^d  therefore  express  D  syinme  trie  ally  in 
terms  of  f   and  r^  ,  namely 

(6.9)   P-^-(V-^r.fr,-.K) 

as  can  be  easily  checked*   Hence  SX    is  the  sum  of  three 
terms, 

(6ao)       Si  ^    K  i-  u  4-  h 


where 


Ok 


V, 

According  to  the  exclusion  principle  ^     must  be 
anti symmetrical  in  the  indices  1  and  2.   Hence  if  we  insist  that 
our  trial  functions  have  the  same  property,  then   ^  J^  will  also 
be  antisymmetrical  so  that  M  -  K.   K  and  L  can  now  be  expressed 
as  surface  integrals  by  Green's  theorem: 


60 


(6.12)        ^ 

where  T,  -  ^\/^\      »      ^»i      is  the  three  dimensional  space   >^2,  -S  ^ 
and   $J:'  is  the  two  dimensional  surface  Aj  ■=.  CI  ,   Thus, 
the  integral  (6,12)  extends  over  a  portion  of  the  five 
dimensional  surface   ^    •   Similarly 

(6.15)        +[    ^^.[^''^r^^f-If  v,^M(^j(d.o 

since,  because  of  the  antisymmetry  of  ^      ,  the  two  parts 
of  L  are  equal. 

We  now  insist  that  a  is  so  large  that  the  following 
two  conditions  are  satisfied:   At  /t.  ^  <^  ,  the  wave  function 
has  split  into  a  function  cp[<?|  \  describing  the  deuteron,  2  and 
3,  and  a  function  'P  [fudescribing  its  motion  relative  to  the 


neutron  I j  i.e. 


(6.14)     ^-  ^l^s)    F(^)  ^^ 


.;\.  >,  a . 


7 
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where  \J      represents  a  negligible  contribution  from  decaying 
modes.   And  further  a  must  be  so  large  that 


(6.15) 


CP(<5^^)      rn       negligible,    ^^  >,    ^  <\  . 


We  restrict  our  trial  functions  to  those  which  satisfy 
corresponding  conditions,  viz. 

(6.14')   "P^    r^     (^4-(^J  ^kM'^^JC     ^     \'>r^ 


and 


(6.15»)     CP^  [^^)     :>     negligible,   f,  ^  |  ^. 

We  assume  that  both   ^(jij  ^nd  ^^\R\^    ^r©  normalized. 
^  I  ^1 )  "^^"^^   have  different  forms  depending  on  which  phase  we 
are  dealing  with.   We  shall  assume  that  its  angular  part  is 
normalized  to  MrTT  .   Thus  for  S  -  scattering, 

^^  ^^  — x^ 

(6.16)  ^ 

and  similarly  for  higher  angular  momenta. 

If  we  neglect  \7  and  V;^  ,  the  integrands  occurring 
on  the  right  hand  sides  of  (6.12)  and  (6.13)  are  expressed 
most  simply  as  functions  of  P   and  Vi  (  see  (6.14)  and 
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(6,14 »)    )    so   that  it   iPi    of  advantage   to  use   these   coordinates 
as   variables    of  Integration.      Therefore   we   transform  from 
IT.     ,    Ty    to    0     ,   Q^    according   to   the   equations 


r,^  IT, 

(6.17) 


which  are  a  consequence  of  the  definitions  (0,3)  and  (0.4). 
It  follows  that 


(6.18)   (/^^(Plf,)  "(t)   (f^^'ljl'^^^  ) 


Also:' 


^OV,  -\%)r,  -^^n) 


(6.19) 

Now  for    /I.-    CL     »     ^      has    the   form   (6.14)    and       <$•  ^     is 

(6.20)  Sf'-  cf[p,)  ^F(n)   -^Sc^(f,)  P(rJ. 

Hence  the  gradient  terms  In  (6.12)  are 

(6.21)  ^t,i  ^  \  cf'(?0  f  (^)?,  -*  ^(?0  f^'(r.)'r, 


63 


(6.22) 


where  f.^Lf^fji^      '  fi  -  f'  /'? i  '  ^""^   '^  Y*^*  )  ;  <JP  I'^J 

denote   the  derivative  of   Ptv"!)  and  J  T^ir,^  along  the  outward 
normal.   Let  us  first  Integrate  over  the  ^^-space,  keeping 
r^   fixed.   The  Integration  extends  over  all  points  with  the 
fixed  r^  and   ^  2^  ^  *   Since  ^yc^  Q.    ,    it  follows  from 
(6.17)  that  at  all  points  outside  this  region  Q.^rr'^       ; 
hence,  in  virtue  of  (6.15)  and  (6.15 ')>  we  may  extend  the 
integral  over  the  entire  P  -space.   Furthermore,  since  CPI^k) 
and   Sc^[p.)   are  independent  of  the  direction  of  f  we  can 
first  carry  out  the  integration  over  the  angle  between  Y^ 
and  O    ,  in  which  all  terms  involving  f^  -  f ,   give  zero, 
while  those  involving  ^^   •  V|   give  MiC  •   Hence 


^-ft)' 


(6.23)      J^Cf^(0  P'(0  ^f*(M  +C^(?0<^1^?^^"^''^^'^J 

(?.)   r^'(^)5F(r,)  -c^UO^'f^?')  P^'^')P'"{n)' 


-'f^ 


^n  'i^.\p''(^)^f'W-^^W<''''(^)];^ 


10 
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due  to  the  normalization  of  Cp(f^).   Notice  that   ^  ^^  uO 
does  not  appear  on  the  surface. 

Let  us  now  specialize  to  the  case  of  S  -  scattering 
and  assume  that  for  ^^^(K    ,      *^:(  Vk)     has  the  form 


(6.24)      F^(0-   A     —j^    -v  « 


where  O^t  is  chosen  to  give  a  smooth  logarithmic  derivative 
across   ^  -^  •  ^  l^\f      ^^^  ^^®  form  (6.16).   Hence 


,JF(a) ..  Jb.  ^ 


u.a 


(6.25)  t\ 

so  that 


^  -  fS 


(6.26)      »-  ^  -.-  »   .  Mil  K. 


Turning  now  to  L  we  observe,  by  (6.19),  that  Vf   acts 
only  on  functions  of  \^\\   and  that 

(6.27)     K/r^  «  [  |twJt^u^  *{  (?,l  )  ^  ?\- 


A    >* 

•  V. 


Hence  the  integrand  contains  ^.  *  m    as  factor  and  the 
integral  over  the  angle  between  9,  and  »i  vanishes. 
Therefore,  by  (6.10)  and  (6.26) 
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Zl 


(6.28)  51^   -i\^     ^-     -   tJTV4   AcJB 


so     that 


(6«29) 


I^^  ^-iri^AB;,    -    ^^^AB 


S  L 


Is  a  stationary  expression  for  tan  h    • 

In  the  derivation  of  (6.29)  v/e  have  met  some  difficulties 
which  did  not  exist  in  the  two  particle  problem:  The  neglect  of 
\?   and  \)    in  (6.14)  and  (6.14*)  and  the  assumptions  (6»15) 
and  (6.15 M  introduce  inaccuracies  for  any  finite  value  of  a, 
even  if  the  interaction  potential  has  a  sharp  cut-off  (  square 
well  );  this  is  in  contrast  to  two  particle  collision,  where 
the  wave  function  has  accurately  the  asymptotic  form  for  all 
distances  larger  than  the  cut-off  distance  of  the  potential, 
and  where  assumptions  of  the  type  (6.15)  ,  (6«15»)  do  not 
enter. 

This  difficulty  can  be  avoided  by  setting  q.  -  oo  but 
then  another  difficulty  is  encountered.   We  have  in  our 
treatment  neglected  all  second  order  terms.   However,  unless 
^  .   has  the  asymptotic  form 

(6.30)     v$^    -^      cf(?.)  f^(r-.)  a^  f^-^^ 
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where   ^|?/  ^^  ^^^  exact  deuteron  function  and  F^  Is 
exactly  sinusoidal  with  wave  number  \<  ,  the  second  order 
terms  are  appreciable  over  an  infinite  portion  of  the 
configuration  space,  and  their  integral,  far  from  being 
negligible,  makes  ^ ^^o   Infinite.   Thus  suppose  that  our  trial 
function  tends  to 

(6.31)     ^;^    — >  I  Lf(?.)  +c}c^(pj^  f>(n)  ,  ^r^^^ 

instead  of  (6.30)  where  T^     is  sinusoidal.  For  large  f.  , 

(6.32)  p     ^    l^r^Jir^'')     -(W:P   -    ^P  J 

where  H  ^  is    the   deuteron  Karailtonian, 

(6.33)  W-p     ^    ^1     X/^^  ^    U^3 

Hence,    from  the  region      ^.    >  -^     ,    say,    in  which    (6.31)    and 
(6.32)   hold  very  nearly,   X^?  gets   the  contribution 


^ 
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where   tp  is  the  expectation  value  of  H  -p  in  the  trial 
state  CPn  +  0  cPp  •   Clearly  AJ  although  a  second  order  term 
(  since   6j>  is  stationary  )  Is  infinite. 

We  are  thus  faced  with  the  following  dilemma: 

(1)  If  we  use  only  an  approximation  for  CO    we  must 
take  a  rather  small  a  or  else  the  effect  (6«34)  which  is 
infinite  for  infinite  a  becomes  important  and  our  neglect 

of  second  order  terms  is  unjustified*   However,  the  finite  a 
introduces  an  inaccuracy. 

(2)  If  we  take  a  infinite  v/e  must  use  the  exact  deuteron 
function  in  the  asymptotic  form  of  ^^   which  may  be  inconvenient 
for  practical  purposes.   (  Such  a  procedure  would  b6  nearly 
impossible  for  more  complicated  collisions.  ) 

The  way  out  of  this  dilemma  is  suggested  by  (6.34),  which 
shows  just  what  the  infinite  contributions  of  the  second  order 
terms  are  and  hence  permits  us  to  subtract  them.   This  is  done 
by  replacing  i<x  ,  (6.6)^X-j  I'^  *  defined  by 


(6.35)   Xol  -.  ^  ^?  ^[^/J(^r,) 


^^ 


whero 


(6.36) 
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Thus  J-ci   depends  in  two  ways  on  the  trial  function: 

(1)  Through  P  whose  form  is  determined  by  the  asymptotic 
deuteron  function. 

(2)  Through  the  trial  -  ^   itself,  when  substituted  in 
(6.35). 

We  shall  show  that 


(6.37)       l'^      +   i2   TT    14    AB^     -     12ti  K  A  S     -t    iS"-^ 


i  ^  ^ 


V 


where      (^      is    the   integral   of  second  order   terms  v;hich  decay 
exponentially   in  all  but   a  finite   region  of   the   configuration 
space. 

Let  us   denote 


'a 

(6.38) 


\^  cf^p'cf  («^r.)Hn)   ^r/fcfl 


Then,   by   (6.36), 


(6.39) 


Now  Vft. 


(6.40)  Xn(_!f)  =   0 
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and,  as  was  shown  above,  without  neglecting  second  order 
terms  we  have 

(6.41)       ^Ifl^C"^)     =   -^  ^>i  Ai^^-B)   ^^ 


where  11m  yU^  -  O 

To  estimate  the  last  two  terms  of  (6.39),  we  decompose 
V^  Into 


where  b  Is  so  large  that  In  V   >      ^     ^J^d  ^  t      have  their 
asymptotic  forms 

(6.43)    ^  =  Cf[^.)f(r,)    ^    $^=  c^i(?.)-r^(rj 


or  the  corresponding  expressions  with  1  replaced  by  2« 
Clearly  the  last  two  terms  in  (6.39)  are  both  of  the  second 
order,  so  that  their  integral  over  the  limited  V^  is  of  the 

^  /I 

type  ^  •  ^   the  other  hand  in  V^  we  shall  find  that  the  ■ 

last  two  terms,  give  zero  owing  to  the  integration  over  fi 

(  or  P^)*      Thus  in  the  region  o^  V   where   H,  >Ji  ,    they  give 
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-Si(*,*-i';^(«,-f;.(6.'-.,)fer](.f,j 

+l^»'-^>')f«^l'>rll»lf,j 


=  0 

Hence 


(6.45)  X,C^^)  -^  J(^p'-^o)|^>('=  A^ 

From  (6.39),  (6.40),  (6.41),  (6.45)  we  now  find 

which  is  the  statement  (6.37). 

This  means  thattin*]  as  calculated  from 

(6.46)  3:'^  4?^  IT  K  A<B^  ::  ^  Tf  R  A'i^V^ 

is  stationary  if  we  admit  only  such  trial  functions  whose 
asymptotic  form  is 

(6.47)    J^  -^  cf^U,)  \  A  ^^:L3t   -t  6^  "ZL^ 
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Thus  the  exact  form  of  the  Internal  wave  function  and  the  • 
binding  energy  of  the  colliding  nuclei  are  not  required. 

Numerical  Calculations. 

We  have  made  exploratory  calculations  on  n-d  scattering 
in  preparation  for  a  more  extensive  investigation.   For 
simplicity  we  have  so  far  only  used  trial  functions  of  group 
structure,  but  we  remark  again  that  this  limitation  is  by  no 
means  inherent  in  our  method. 

V/e  have  assumed  an  ordinary  Wigner  potential 

(6.48)    V.I)  •=-  Vj,  <2:     "^^ 

with  Vo-"  72.00  mc   and  a  ■=.  0.8  in  units  of  A^,  c^   2.80  x  10  cm. 
These  constants  give  the  correct  binding  energy  of  2.18  MEV. 
for  the  deuteron. 

As  trial  functions  for  S-scattering  we  have  used 

(6.49)     ^^  =  Cf^(?J^H(«'0'X, -^c^^(?0-f^('r07^   ^ 

where  we  have  set 

43«   These  constants  are  almost  identical  with  those  of  Motz  and 
Schwlngor,  ref.23,  who  however  used  an  exchange  potential* 
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(6.50  cont»d) 

'Y I  -  -  T(  j^  r:  (J^[\)    d[l)    o^[^)  dw^OAXei  Sc^cMMa.^ 

All  distances  are  measured  In  units  of  Aq    •   ^tT  ( ?i)  ^^  ^^® 
best  approximation  of  its  type  for  the  deuteron  function. 

The  polynomial  approximation  for  ^l*^)  has  the  v/rong 
asymptotic  form.  Hence  it  was  assumed  to  hold  only  inside  the 
region  A  <  C\.  . 

We  have  calculated  only  S-scattering  at  zero  energy.   Hence, 
asymptotically 

(6.50-)    -^  [n    -^    f^  -i  ^ 


where  A  and  B  are  chosen  so  that  the  polynomial  form  of  ■V,^[«'j 
passes  smoothly  into  (6.50*)  a.t     j\ -z^  (\.   • 

As  we  are  limiting  ourselves  to  a  finite  region,  the 
difference  between  l(a)  and  I* (a)  is  negligible  (  see  (6.55) 
and  (6.36)  ).   The  calculations  to  be  described  were  carried 
out  by  using  I  Ma)  throughout. 

We  have  calculated  limit  (  tan  ^/  \^  )  ~  X,  for  quartet 
and  doublet  scattering,  which  determine  the  total  cross  section 
in  the  limit  of  zero  energy  by  the  equation 
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(6.53)     6*  ^  4u  (^  3  X«\   -^  -^  ?^*l  J 

where  q  and  d  refer  to  quartet  and  doublet  scattering 
respectively. 

In  the  case  of  doublet  scattering,  the  convergence 
with  increasing  values  of  v\  was  slow.   Since  in  the  case  of 
two  particle  scattering  method  2  gave  somewhat  better  results 
than  method  3,  (  see  Table  I,  p.  30  )  we  have  calculated  X  j 
not  only  by  (6.29)  but  also  by  the  equation 

8  ^ 

which  bears  the  same  relation  to  (6.29)  as  does  (2.32)  to 
(2.23).  X  is  determined  from  the  logarithmic  derivative  L 
according  to  the  equation 


^L 


\/      ^ 

(6.52)     X.  ts  -  - 


Hov/ever,  the  convergence  was  only  slightly  improved. 

The  approximate  values  of  X,  obtained  by  our  variational 
methods  depend  both  on  a  and  n.   For  large  enough  n  we  expect 
the  following  variation  with  a: 

(Al)  When  a  is  too  small  the  expression  (6.29)  for  tan 
is  incorrect.   Hence  X  should  depend  strongly  on  a* 

(A2)  When  a  is  sufficiently  large  for  (6.29)  to  hold, 
X  should  be  fairly  independent  of  a  until 
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(A3)   a  becomes  so  large  that  P^  can  no  longer 
approximate  the  true  function  over  the  entire  interval.   Our 
variational  principle  than  breaks  down  because  second  order 
terms  become  important*   Here  again  we  expect  a  considerable 
variation  of  X  with  a. 

Thus  we  expect  to  find  an  extremum  in  the  region  (A2) 
and  the  value  of  X  at  that  extremum  may  be  taken  as  our 
approximation  for  X. 

As  we  increase  n  we  expect 

(Bl)   Gradual  extension  of  the  "plateau"  (A2)  towards 
larger  values  of  a. 

(B2)   Increasingly  smaller  differences  X  (  a;  n  +  1  ) 
-  X  (  a;  n  )  for  large  enough,  fixed,  a. 

In  an  approximation  method  n  should  be  increased  until 
(Bl)  and  (b2)  are  so  well  satisfied  that  further  increase  of 
n  is  not  likely  to  change  X  appreciably. 

In  virtue  of  the  choice  of  our  potential  and  trial 

function,  all  integrals  occurring  in  I^^  could  be  analytically 

are 

evaluated.   Our  results  for  Xa  and  X  j^given  in  Tables  III, 

IVa  and  IVb  and  plotted  in  Graph  2.   The  best  3  -  parameter 
wave  functions  r^  fr)  determined  in  the  usual  way  from  our 
variational  principles  are  shown  in  Graph  3. 


GRAPH  2.  NEUTRON  g^UTERON  SCATTERING 
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Table    III. 


Quartet  Scattering   ;   X  g  (a;n)   /Aq  , 
(    Method  3,     (6«29)    ) 


0 

1 

2 

2,5   . 

-    4.39 

--    2.16 

-   1.58 

3,0 

-   4.39 

-    1.91 

-   1.85 

4.0 

-  4.39 

-    1.84 

-    1.83 

6.0 

-   4.39 

-    2.03 

-  1.84 

8.0 

-   4.39 

-    2.26 

-   1.92 

Table  IV  a. 


Doublet  Scattering  ;  Xj^(a;n)  /^o 


(  Method  3,  (6.29)  ) 


^^^v^l 

0 

1 

2 

a/Aa\ 

2.0 

+   8.53 

- 

1.18 

- 

1.58 

2.5 

+   8.63 

- 

2.03 

- 

3.73 

2.6 

+   8.53 

- 

2.14 

- 

1.03 

3.0 

+   8.53 

- 

2.51 

- 

1.58 

4.0 

+   8.53 

- 

3.38 

- 

2.01 

6.0 

+   8.53 

- 

5.65 

- 

2.98 

8.0 

+   8.53 

- 

9.16 

- 

4.05 
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Table  IV  b. 


Doublet  Scattering  ;  X^  (a;n)  /  ^o  • 


(  Method  2,  (6.51),  (6.62)  ) 


^^^J\ 

0 

1 

2 

a/^S^^-^ 

2.5 

-   3.53 

-   0.36 

-    1.51 

3.0 

-   4.62 

-    2.04 

-    1.47 

4.0 

-   7.53 

-    2.85 

-    1.90 

6.0 

-   20.23 

-   4.47 

-    2.78 

Interpretation  of  the  numerical  results: 

A.  The  convergence  of  X^  is  much  more  rapid  than  that 
of  X  J  .   This  is  due  to  the  fact  that  P(  r*  )  is  more  irregular 
for  doublet  scattering  (  see  Graph  3  )  and  therefore  requires 
more  terms  for  adequate  approximation. 

B.  The  variation  of  X  a  with  a  and  n  is  as  expected 
(  pp.  73  and  74  );  X  ^^  seems  to  behave  in  a  similar  manner 
although  the  number  of  parameters  is  not  large  enough  to 
exhibit  a  proper  plateau. 

C.  To  determine  the  total  cross  section  we  must  select 
the  best  values  of  X  <3  and  X  Jl  from  our  tables.   In  the  case 
of  X  a  (  Table  III  )  we  have  evidently  come  very  close  to  the 
limit  which  can  be  reached  with  group  structure  form.   We  see 
that  our  results  depend  markedly  on  a  for  values  smaller  than  3* 
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This  indicates  that  the  interaction  is  appreciable  out  to 

about  that  distance.   We  therefore  choose  X  a  =-  X  ^.  (4;2)  -  -1«85^q, 

The  choice  of  X  ^  is  more  difficult.   At  a  =  3  we  cannot  be 

sure  that  the  interaction  is  completely  negligible;  on  the 

other  hand  at  a  r.  4  we  are  still  quite  far  from  convergence. 

Let  us  nevertheless  choose  X^-  X  ,(4;2)  -^    -  1.90  a^    (  see 

Table  IV  b  ).   This  gives,  by  (6.53),  ^ --.     3.39  barns.   If 

we  estimate  the  limit  of  X  jl  which  could  be  reached  by  taking 

larger  n,  as  X  ^  t^   -  l.SOyi^we  obtain  ^  =-  2.94  barns.   These 

figures  are  in  qualitative  agreement  with  the  experimental 

value  of  3.2  barns 

Let  us  summarize  the  possible  inadequacies  of  this 
tentative  calculation. 

(a)  Wigner  forces  may  not  be  adequate,  to  describe  this 
process. 

(b)  Our  approximation  for  Cf  [?i\  (  s^®  (6.50)  )  may 
not  be  sufficiently  accurate. 

(c)  The  group  structure  approximation  may  not  be 
adequate. 

(d)  The  Schroedinger  theory  may  not  be  capable  of  an 
accurate  description  of  this  process. 

Of  these  shortcomings  (a),  (b)  and  (c)  can  be  remedied 
by  more  extensive  calculation  of  the  same  kind.   When  this  is 
done,  it  will  be  possible  to  decide  to  what  extent  (d)  is  true» 

In  future  calculations  the  following  points  should  be      ' 
given  consideration! 

(1)   Terms  allowing  for  the  polarization  of  the  deuteron 
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should  be  introduced. 

(2)   It  may  be  of  advantage  to  use  trial  -  P's  of  the 
correct  asymptotic  form.   The  reason  for  this  is  the  following: 
If  we  use  a  polynomial  expression  for  F  {^) ,    we  must  include 
a  large  number  of  terms  since,  due  to  the  loose  structure  of 
the  deuteron,  the  region  of  interaction  is  very  large 
{'^C       8  X  10   cm  ) .   This  necessitates  the  calculation  of  a 
great  many  matrix  elements,  which  is  probably  more  laborious 
than  the  numerical  integrations  which  become  necessary  v/hen 
functions  of  the  correct  asymptotic  form  are  used. 


!>?♦   The  General  Case  of  Composite  Collisions. 

The  method  used  in  the  preceeding  section  can  be 
extended  to  all  collisions  of  composite  nuclei  in  which  triple 
or  multiple  disintegrations  are  energetically  forbidden.   The 
analytic  aspects  of  the  problem  are  quite  analogous  to  those 
encountered  in  n-d  scattering.   However  the  combinatorial 
side  of  the  general  problem  is  somewhat  complicated  and  no 
detailed  description  of  the  general  case  will  be  given  here. 
We  shall  merely  mention  a  few  points  of  interest. 

Modes  of  Collision. 

Our  picture  is  the  following:  We  have  a  number  of  pairs 
of  nuclei  colliding  and  a  number  of  pairs  (  possibly  different 
ones  )  emergin  •  after  the  collision.   The  groupings  of  nucleons 
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Into  pairs  are  alternative  (not  independent)  processes, 
each  pair  comprising  the  same  nucleons.  There  is,  then  a 
"region  of  interaction"  in  the  configuration  space,  where  all 
nucleons  are  close  together,  and  an  "asymptotic  region", 
where  the  nucleons  have  separated  into  pairs. 

On  a  hyper-surface,  2^  ,  separating  the  two  regions, 
the  wave  function  has  the   form 


(7 


.1)  4r  -  1  $ 


(I )      iv ) 

where,  if  (O.    and  cP;  describe  the  normalized  internal  wave 

functions  of  two  colliding  nuclei  and  P  their  relative  motion. 


(7.2)   $.  :.   C|>,  C^,   P"-. 


A  function  of  the  type  (7,2)  will  be  called  a  mode  of 
collision. 

If  S^  is  a  proper  function  of  a  complete  set  of  operators 
(  H,  J,  P,  JItj..),  its  modes  satisfy  the  following 
asymptotic  orthogonality  conditions: 

(7.3)  \  ^;<5^  1^  *0     ^  ■^^■^. 
t 

44.      The   terra  "collide"  will  be   used  in  the   remainder  of  this 
section   to  describe   both  collisions   and  disintegrations. 
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For  if  Op  t    and  Of;  describe  different  pairs  of  nuclei, 
they  exist  on  different  portions  of  the  hypersurface  Z.  and 
hence  (7»3)  is  satisfied.   If  however  they  describe  the  same 
pair  of  nuclei,  then  either  cP-   •*  Cf  j  or  <^;  4^),  so  that  the 
orthogonality  (7«5)  follows  from  the  integration  over  the 
internal  coordinates  of  the  nuclei. 


Number  of  Arbitrary  Constants  of  the  Scattering  Matrix. 

We  have  remarked  in  section  1.5.  that  due  to  its  unitary 
character  and  symmetry,  the  two  dimensional  scattering  matrix 
involves  only  three  independent  constants.   In  the  present 
case  the  scattering  matrix  has  the  same  order,  n,  as  there  are 
different  modes  of  collision.  Wcien   represented  in  the  scheme 
of  incoming  and  outgoing  modes,  we  have  proved  that  the  number 
of  independent  constants  is  (n)  (n+1) .   Only  this  many  constants 
need  to  be  determined. 

r 

Rayleigh-Ritz  Method. 

7/e  have  already  applied  the  Rayleigh-Ritz  method  to  two 
problems  involving  two  modes  of  collision:  n-p  scattering  v/ith 
tensor  forces  (  section  1.5.  )  and  n-d  scattering  (  section  1.6. )• 
In  section  1.5.  we  have  discussed  two  different  variational 
methods:  One  dealing  with  complex  (  running  )  wave  functions 
and  yielding  all  elements  of  the  scattering  matrix;  the  other 
dealing  with  standing  waves  and  giving  only  the  proper  values 
of  the  scattering  matrix.   (  Only  this  second  method  was  used 
In  n-d  scattering.  ) 


81 


These  methods  can  be  generalized  as  follows: 

Running  waves : 

When  dealing  with  a  special  case  of  running  waves  in 
section  1. 5.,  we  introduced  a  wave  function,   vl^^   , 
(5«5),  whose  angular  and  spin  functions  were  the  complex 
conjugates  of  those  of  the  function   vK  •   This  was  done 
in  order  to  make  use  of  the  orthogonality  relation, 
corresponding  to  (7«3) 

In  the  present  more  general  case,  such  a  procedure 
is  not  convenient,  since  the  v/ave  function  does  not  everywhere 
separate  into  radial  and  angular  parts.   On  the  other  hand, 
if  v/e  take  the  complex  conjugate  of  the  entire  function,  the 
sense  of  the  incoming  and  outgoing  waves  is  reversed  and  the 
(pt  *s  and  n  ^s,  denoting  their  amplitudes  are  mixed  up.   To 
avoid  this  difficulty  we  use  the  following  notation: 

(7.4)     ^^-r^  ^   V^l>^)^ 

where   J    and  *$"  '   have  the  asymptotic  forms 


(7.5) 


S  •  is  the  ane;ular  function  of  relative  motion  of  the  mode 
We  note  that  by  the  definition  (7«4), 


^-^^    -\M)^  ^-^)     ^)^ 


(7.6)  at'        :.  ft: 
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which  Is  the  reason  for  introducing  i   ,  as  only  in  this  way- 
can  the  oL   's  and   p  ^s  remain  associated  with  incoming 
and  outgoing  waves  respectively. 
We  now  form  the  expression 

where   x   and  ^   have  the  same  quantum  numbers,  V  is  the 
region  of  interaction,  H  the  total  Hamiltonian  and  E  the 
total  energy  of  the  system.  The  first  variation  of  I^,* 
can  be  transformed  into  a  surface  integral  in  the  usual 
way  and  we   are  led  to 


(7.8) 


^    AU/i.  >w<X^S 


(Compare  (5.7))  Therefore  the  following  equation  is 

C7>t*»    .net  J 
stationary  relative  to  variations  of  3^  and   zS 


We  would  like   to  convert   this  into  a  stationary 
expression  for  the  elements  of  the   scattering  matrix 
(see  (5.9)).  Now  in  the  case  where  groups  may  collide 
with  different  kinetic  energies,  we  must  modify  the 
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definition  of  the  scattering  matrix  given  in  (5.1).  The 
amplitudes   A'  are  linear  combinations  of  the  OCj  ,  viz. 


(7.1U) 


f- 


f';    --?  '-)■'■,  ■ 


The  total  flux  into  the  region  of  the  conf igiiratlon  space 
enclosed  by  ^  is 

A  * 


(7.11) 


m4 


)         l*) 


The      last   line   shows  tiiat    the  matrix 


(V.12) 


is  unitary,  and  hence  may  properly  be  taken  to  be  the 
scattering  matrix.  If  in  analogy  with  (5.2)  we  define 


I  r= 


(7.13)     (^  c 


K,  o<, 


^t  «<t 


%'. 


\ 


46.      A»   Wlntner,    Spektraltheorie   der  Unendllchen  Matrlzen, 
(Leipzig,    1929)    p.    34. 
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then  by.^7,lu)    and    (7.12;, 
(7.14)  1?>    m     S  •   ft  . 

In   this   notation    i7,9)    can  be   written     as 

,7.15;   X,,  ^  -^  A   .  B  ^.^  "  T  ^   ^ 

Hence  by  choosing  our  trial  functions  with  only  a  single 
non-vanishing  component  of  A,  we  can  obtain  a  stationary 
expression  for   any  one  matrix  element  of  S  at  a  time. 

Standing  Y/aves: 

The  modes  of  the  proper  functions  of  the  scattering 
matrix  all  have  the  same   phase  shift  (see  section  1»5,  ,(^.\^)'^^.) 
We  form 


X  ^  $  vf  ^[^u  +e)  Sk'    alt: 


(7.16) 

_       V 

where  3^  is  a  trial  function  v/hich  consists,  asymptotically, 

of  standing  v/aves.  The  first  variation  of  X  now  becomes 

(7.17)    sr---il^:.  ^;  <^^t^1; , 

where   j^  •  is  the  amplitude  of  %y^  ^j^A^  jys'     in  F  ,  the 

function  of  relative  motion. 

Hence 
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(7.18) 

The  collision  determinant,  which  arises  on  application 
of  the  Rayleigh-Ritz  method  is  of  the  Y\       degree  corresponding 
to  the  sum  of  n  squares  multiplying  tan  y\   .   Its  n  solutions 
are  approximations  for  the  tangents  of  the  n  proper  values  of  . 
the  phase  shift. 
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II.   THE  VARIATION  ITERATION  METHOD  IN  COMPOSITE  COLLISIONS. 


II. 8.   Preliminary  Considerations. 

Schwlnger   has  recently  developed  a  variation  iteration 
method  for  two  particle  collisions.   His  treatment  leads  to  a 
stationary  expression  for  tan  n  ,  involving  a  Green's  function 
and  giving  very  good  results,  even  with  rather  crude  trial 
functions.   In  this  section  we  shall  first  rederive  his 
results  in  a  slightly  different  way,  to  show  its  relation  to 
the  bilinear  variational  principles  discussed  in  Part  I;  and 
then  shall  examine  this  method  with  a  view  of  generalizing  it 
to  many  particle  problems. 

Consider  first  the  one  dimensional  Schroedinger  equation. 


(8 


•"  l£.-'^'*'>-(wi-«i-^ 


where   u  satisfies    the   boundary  conditions 

(8.2)       Lc[o)aO-     u,[x)     ^    SU  id^  -vtax^ii,   c^,  K/^ 

Let   the   Green's  function  .G(x,   x')    satisfy  the   differential 
equation 
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and  certain  boundary  conditions    to  be    determined  later. 
Multiplying    (8,1)    by  G    and    (8.3)    by  u,    and   integrating   over 
X   from  0   to   a  we   find 

(8.4)  ^ 


To  make  the  integrated  terms  vanish  vie   now  impose  the  same 
boundary  conditions  on  G  as  ar^  satisfied  by  u,  namely 

(8,5) 

where   F    (xM    is   some   function  of  x*.      The    solution   of    (8#3) 
and   (8.5)    is 

(8.6)  <^  V^;  ^/     ~     —     Sv^  K.  3[:^  .  6ik5  KX;,  ^     ~--4-  S w  K. :? .  ^ C^ Ko: 

1^  Is 

In  virtue  of  (8.4)  and  (8.5)  we  now  obtain,  on  letting  Cc-^oO, 

(8.7)  i^[:c)  -  1    \    GU:r')  ^(x')  uU')  ^'. 

0 

This  may  be  regarded  as  a  proper  value  problem  for  A-^ 
for  a  given  cotw  .  and  the  discussion  can  be  completed  from 
this  point  of  view:   (8.7)  Is  multiplied  by  f (x)u(x)  and 
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integrated   over  x,    giving 


0 

(8.8) 


o^ 


+7^{Uj<aW  c„uv)  {[v'j<^u;j  ju  <u 


where 


(8.9)  G^\jii^y})     rz,      J^   Si*^^  k^^.  ^^  ^^>  . 

(8.7)  is  a  stationary  expression  for  7-  corresponding  to 
a  given  cot  wi  or  for  cot  Wi  corresponding  to  a  given  A  • 
On  the    other  hand  we  may  derive  (8.8)  as  follows; 
We  note  that  u  defined  by 

(8.10)  W^^%\(^   U.^'j  ^[^'}  l^Cx')  oU'^ 


becomes,  asymptotically 

(8.11)  ^ 


li::^f  s;^icx'^(v)u(.'i^/ 


K   * 

0 

and  thus  has  the   correct  asymptotic  form  and  phase  shift 
for  any  trial  u;  furthermore  by  (8.7),  u  j^x  u  when  u  is 
the  correct  solution.  Hence  if  used  in  our  bilinear 
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variational  principle    (2,48), in  conjunction  with  u    =.    u, 
we  have    simply 

.  O  ; 

for  arbitrary  variations  of  u.  This  is  of  course 
equivalent  to  (8.8). 

Three  Dimensional  Formulation^ 

Turning  now  to  the  solution  of  the  equation 

(8.15)  [v'-t  k'^-V  >-((0]4'(''^  =  '^- 

one  can  proceed  in  an  entirely  analogous  manner*  The  solution 
of  (8.13)  corresponding  to  a  wave  of  unit  amplitude  incident 
along  |4   is  written  in  the  form 

and  becomes  asymptotically 

Hence  the  scattering  amplitude  in  the  direction  K.  is 
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given  by 

(8.16)   ^(k^,k.)  -.'~\^''''^\i^')^,A^')[^r') 

This  permits  us  to  write  (8.14)  in  the  homogeneous  form 


(B.17)   ^^J,J  ._  ';\  j(r(r.r')  {{r')  ^^,ir')    (^'') 

where 

For  large  values  of  r,  (8.17)  becomes,  in  the  direction  |< 
(3.19) 


{i^'J4'xao(<i^') 


Hence  no  matter  what  function  v>   /  f^\   is  substituted 

in  the  right  hand  side  of  (8.17),  the  left  hand  side  always 

has  the  correct  scattering  amplitude  f(  k^i.j^|)» 

This,  in  fact,  is  the  reason  for  replacing  (8.14)  by  (8.17). 

Returning  to  the  variational  principle,  it  is  clear 
that  if  we  use  only  such  ^^   which  have  the  form 
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and  arbitrary  ^^  ^r)  ,  then  since   -P;f  L^hj  ^x)    =  -([v^l^kJ 
it  reduces  to  ^  ^^    u  -0*  V/hen  this  is   written  out  in 

full,  it  serves  as  a  variational  principle  for  the 
determination  of  -{(Kj^.kJs 


(8.21)  - 


t-K,lO-((0  C^kJO^^k) 


Let  us  note  specially  that  no  boundary  conditions  had  to 
be  imposed  on  either  ySCu,    or  ^K,  •  This  is  in  keeping 
with  the  fact  that  all  integrals  occuring  in  (8,21)  extend 
only  over  the  region  of  interaction,  where  f(r)  and  f(r*) 
are  not  zero. 

The  one  and  three  dimensional  treatments  have  in 
common  certain  features  which  may  be  expected  to  carry 
over  into  a  similar  treatment  of  composite  scattering 
processes: 

A*   Both  treatments  involve  one  iterationnS.lO) 
and  (8, 20 )j  and  hence  give  excellent  results  even  with 
poor  trial  functions.  In  principle  the  accuracy  can  be 
increased  by  iterating  a  number  of  times. 
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B,   The  trial  functions  need  not  satisfy  any  boundary 
conditions;  instead, 

G.   The  boundary  conditions  are  automatically  satisfied 
by  a  proper  choice  of  the  Green's  function. 

D.   The  Green's  function 

!•   except  for  the  inhomogeneous  term,  satisfies 
the  same  differential  equation  as  does  the  wave  function 
asymptotically  (  see  (8«3)  and  (8.18)  )  and 

2,   has  the  same  phase  shift  (  or  scattering 
amplitude  )  as  the  correct  wave  function. 


II.9.   Neutron  Deuteron  Scattering;   Six  Dimensional  Green's 
Functions. 

As  a  prototype  of  composite  scattering  processes  we 
shall  consider  in  this  section  the  collision  of  a  neutron 
with  a  deuteron.   We  shall  generalize  the  method  described 
in  the  last  section  in  a  straightforward  manner,  but  shall 
find  that  for  composite  scattering  it  is  too  complicated  to 
be  of  practical  use. 

The  wave  equation  for  our  system  is,  by  (6.3)  -  (6.5) 

(9.1)     T^E    (-C^-u^^-o^^^a^^-tK^eoj^-^O^ 


where 
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(9.2)  D--Cv-.^v/^:)--(vr:^ivj:). 

Let  us  for  the  moment  ignore  the  spin  dependence  and 
simply  insist  that  the  wave  function  be  anti symmetrical 
in  1  and  2.  Then,  asymptotically  the  S-wave  has  the  form: 

(9.3)  ' 


^v 


The  wave  equation  (9.1)  tecomes  asymptotically: 


(9.4) 


According  to  D»l«,  section  II. 8,  the  Greenes  function 
should  satisfy  the  differential  equation 

(9.5)     L  ^   ^  -  S(r»-r,')  S  [h-fi] 

in  the  asymptotic  region  of   the  configuration  space, 
where  L  reduces  to  L  and  L  for  large  values  of  n.  and  ^i. 

\        2.  I 

respectively.  In  the  one  dimensional  problem  L  was  simply 
3-  -i-  K     ^^^   made  no  reference  to  the  interaction 
potential;  hence  a  solution  for  G  could  be  easily  found. 
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The  fundemental  difficulty  In  the  case  of  composite 
collisions  is  that  parts  of  the  interaction,  namely  the 
internal  forces  of  the  colliding  nuclei,  persist  at  infinity^ 
Therefore  they  appear  in  the  asymptotic  operators  L  and  L 
and  hence  also  in  L,  v/hich  makes  the  determination  of  G 
practically  impossible.  In  fact  the  simplest  L  which  we  can 
choose  is  L  =:  P  +  "^n.*  since  v  ^^^  is  the  only  interaction 
which  disappears  asymptotically.  Of  course  the  practical 
solution  of 

(9.6)     (^_D  -u^^-  l/jj  +«"-+  ^t)  C  -  ^^iu-fi')  ^("i.-f^J 
is  quite  impossible. 

For  this  reason  v/e  shall  not  pursue  this  method  any 
further,  but  merely  state  that  If  formally  carried  through, 
it  has  the  properties  A,  B,  C  and  D  mentioned  in  the  last 
section  and  leads  to  stationary  expressions  for  the  elements 


of  the  scattering  matrix  (  1  and  J  refer  to  the  incoming 
and  outgoing  neutrons,  V4  and  Ka  to  their' directions)  • 

Another  formulation  v/hich  suggests  itself  is  analogous 
to  the  treatment  of  the  bound  state:  We  write  the  wave 
equation  as 


and  introduce  a  Greenes  function  to  satisfy 
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(9.9)  (yrM^?:j>C(r3,?3  |<,fO=-<^(^-*-3'j^(f.-j3'j. 


A  solution  of  this  equation  is 

ri'  Hit 


"•"'  '^Mrr;^-t4^°t^"-'3J^^9 


where 


(9.11) 


Since   ^  ^  Q-p   is  negative  in  our  case,  G  falls  off  exponentially, 
Hence,  on  multiplying  (9,8)  by  G,  (9,9)  by  4^  ,  subtracting 
and  integrating  over  r     and   P   ,  the  surface  integral 


vanishes  and  we  are  left  with 


However  the  exponential  vanishing  of  G  wipes  out  all  reference 
to  the  boundary,  from  where  in  the  past  we  had  always  derived 
the  connection  with  the  scattering  properties.   Thus  the 
present  formulation,  while  providing  an  iteration  procedure 
for  the  collision  wave  function,  does  not  lead  to  a  stationary 
expression  for  the  phase  shift  or  scattering  amplitude. 
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11.10.   Neutron  Deuteron  Scattering;   Three  DlTnenslonal 
Green's  Function. 

'  Since,  as  we  have  seen^the  straightforward  generalization 
of  the  two  particle  procedure  is  of  no  practical  use  we  have 
tried  the  following  alternative  approach: 

Where  fj-^-ao,  the  wave  function  separates  and  therefore 
satisfies  not  only.  L.St^^^    ^  ^^^  (9.4)  )  but  also 

(10.1)  r,  -^  ^  •C^r,'-^^^^)  '^  C?.>  O  -0 

Hence,  in  partial  agreement  with  D.2  (  section  II. 1.  )  we 
introduce  a  three  dimensional  Green's  function,  which 
satisfies 

(10.2)  ^r.^'-tK'^)  ^(rijr/)  r  _  S(vf;) 

and  has  the  correct  scattering  amplitude  in  the  direction 
y^^  •   Suppose  that  we  consider  the  wave  function,  L^  yJi^ 

corresponding  to  a  plane  wave  of  neutron  1,  of  unit  amplitude, 

ft 

incident  along  [i^   while  both  neutrons  1  and  2  are  scattered  • 


Then  the  scattering  amplitude  of  neutron  1  along  Ka    is 
Of^t'J   K^^M  •   A  Green's  function  satisfying  (10.2)  is 

46.   This  is  in  conflict  with  the  exclusion  principle  but 
convenient  for  our  present  considerations.   By  super- 
posing the  symmetrical  wave  function,  with  a  negative 
sign,  a  solution  obeying  the  exclusion  principle  can, 
of  course,  be  obtained. 
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(10.3)  Hu  [u-Y^'X 


.     u(*i     .  .      .^^)  » 


It   satisfies    the    differential   equation   and  asymptotic 
conditions   of   the   wave   function  only  where    ^\  — ^  o^    ,    but 
not  where    T^  — 9>  cO  ,    and   thus    not   over    the    entire    asymptotic 
space    (    see   section  II. 1,,    D.l.    )• 

Proceeding  nonetheless  we   v/rite   the   wave   equation    (6.5) 
for   our   system  in   the    form 

where 

(10.5)      -U.(.<?o^)   H    -^  \7^%   U^jH-Uj,    +Ua   -6o 

is  a  kind  of  effective  potential,  acting  on  particle  1  when 
f*^  is  small. 

From  (10.2),  (10.3)  and  (10.4)  we  can  deduce  in  the 
usual  way  that  the  wave  function  satisfies  the  integral 
equation 

(10.6)  ^Kii'buM  -^(^(n;v.')au.,^'j^,7(?.,r:j(5ir;] 

as  may  also  bo  checked  directly.   For  large  r^  in  the 
direction  ^^    and  small  values  of  ^.    this  becomes 
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^i 


^^^'1%'d  -^  ^ — T""     [^       ■+ -f(^w,«*3^- 
(10.7) 

Thus  again,  regardless  of  the  \l^^'J(^j  r,M  used  in  the  right 
hand  side  of  (10.6),  the  left  hand  side  has  the  correct 
scattering  amplitude  along     «   However  unless 
splits  for  large  )\     into  a  product  of  the  form 

(10.8)    ^c(?o^')  -^^(f.)  ?^(r;j 

LL  j^u^*^    does  not  vanish  so  that  (10.6)  would  diverge. 
For  in  contrast  to   ^ilviULis  not  a  simple  multiplicative 
factor  and  "cuts  off"  only  wave  functions  of  the  type  (10.8). 
Now  consider  the  expression 

(10.9)  I-,  \    'i.^'>[<^^^^ii''-^^[^„r.))'^^^''i'^U)W 

Clearly,  for  the  correct  *^Kj  ,1^0.  Let  us  now  consider 
the  variation  of  I,  if  we  admit  only  trial  functions  of  the 
form 


(10 


•  10)    ^K^^  -^  \  G  tsO  U  U; '')  5c(f'>^'X^^0 


where       \B\<^'^will  be   suitably  res  trie  tedj    and      ^^^f*'   is 


unrestricted.      Then 
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(10.11)    SI  -^i^^-«ri-'0)s'5',--  S^,>.;(-p)^-K^'}t<r.%) 

As  in  section  1.6. ,  (6.8)  -  (6.13)  this  leads  to  surface 
integrals  over  the  five  dimensional  surfaces  ^  ^^       [  '^ i  -  ^  ) 
and  2- a  l-'^i-^)'  ^^®  *^^®  over  2L  ^  vanishes  because  ^^ 
has  the  correct  Ai^oL  )  ^(^    1  owing  to  its  expression  in  ter 


of  a  Green's  function.   The  one  over  Z«  vanishes  if  S  ^  ^ij^^) 
describes  only  outgoing  waves  of  neutron  2,  since   x-^<f^ 
contains  only  such  waves  and  hence  the  integrand  vanishes. 

We  shall  show  that  ^ti^^')  has  indeed  this  property 
(  of  outgoing  neutron-2  waves  )  if  only  fc^-*  has  it.   For 
^j^zH      ^^^    ^"^i^Qf  ?    exceeds  2/3  a. 
Hence,  by  our  assumption. 


(10.12)   <$*  K^'  -~    c^(fi)-  -f 


Now  for   Q  >y     2/3  a,    (10.12)    is   an  exact  solution  of   the  wave 
equation   (10.4)    so   that 


(10.13)     [v,;-tK"^   U(f.O]   ^ 


xr   =  ^ 


Hence 


(10.14) 


p  im>:'^tj*rm'mmM$ 
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(10.14  cont'd) 


•^,'=A 


But  when  P  ^  2/3  a  and  T,  =  t\  ,   ^  kJJ'  (,  fu  "^I  J  ~  ^  ' 

so  that  the  surface   Integral  vanishes.  Hence,  using  (10.8), 
we  find 

as  was  to  be  shown. 

Therefore  if  the  trial  function  ^  ^^]>     satisfies  the 
conditions 

(10.16)  <p  ^0)  ^ 


the  equation 


0.1  =  [*.„«(,,,,)  a(5,r,J$,„  (f„r,)(^J(J,J 
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provides  a  stationary  expression  for 


iWjju'P- 


In  this  way  we  can  obtain  the  amplitude  for  elastic 
scattering  in  all  directions.   In  the  simple  case  of  n  -  d 
scattering,  the  amplitude  for  exchange  scattering  is  then 
determined  owing  to  the  restrictions  to  which  the  scattering 
matrix  is  subject.   Thus,  let  us  consider  a  sub  -  matrix  (^^;  J 
corresponding  to  scattering  with  a  given  angular  momentum,  L 
where  the  subscripts  refer  to  the  incoming  and  outgoing 

neutrons.   There  are  two  modes  of  disintegration  so  that,  by 

-  u 

section  1.7.  (  j^'  )  involves  at  most  three  arbitrary  constants. 

In  fact,  due  to  the  equality  of  the  two  neutrons  this  number 
is  further  reduced  to  two.   Thus  a  single  diagonal  element  com- 
pletely defines  the  matrix.   It  is  clear  therefore  that  the  • 
determination  of  ^(k^^I^aI,  from  which  all  elements  aj(  ,  can 

be  deduced  is  sufficient. 

A  (Ail     iv\ 
Variational  principles  for  -Ki^^^i }  ^^ )    can  be  obtained 

in  a  similar  way,  but  in  addition  to  volume  integrals,  as 

occur  in  (10.17)  they  also  involve  surface  integrals.   The 

same  holds  for  variational  expressions  for  the  proper  phases 

of  the  scattering  matrix. 

No  numerical  application  of  (10.17)  or  its  equivalent 

for  the  phase  shifts  has  been  made.   The  labour  involved 

would  be  enormous.   Thus  the  calculation  for  the  phase 

shifts  involves  a  six  dimensional  integral  which  can  be 

reduced  to  five  dimensions  by  an  analytical  integration. 

The  remaining  integral,  owing  to  the  presence  of  the  Green's 

function  must  be  evaluated  numerically.   Furthermore,  since 
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the  iteration  (10.10)  extends  only  over  three  dimensions,  it 
is  not  certain  how  much  it  improves  the  trial  function  and 
hence  one  cannot  be  sure  that  (10.17)  will  yield  very  good 
results. 
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CONCLUSION. 

The  problem  of  many  particle  nuclear  collisions  is  a 
difficult  one.   Even  in  classical  mechanics  the  treatment 
of  several  interacting  particles  is  in  general  complicated 
and  hence  one  cannot  reasonably  expect  a  royal  road  to  the 
solution  of  the  corresponding  nuclear  problems. 

As  remarked  earlier,  however,  good  results  have  been 
obtained  for  many  particle  bound  states,  with  a  reasonable 
expenditure  of  effort  and  time.   In  Parts  I  and  II  we  have 
examined  the  tv/o  methods  which  have  been  used  for  bound 
systems,  for  their  applicability  to  collision  processes. 

We  believe  that  the  Rayleigh-Ritz  method,  discussed 
in  Part  i,  provides  a  practical  procedure  also  for  the 
solution  of  collision  problems.   To  be  sure,  the  fact  that 
the  trial  functions  must  satisfy  certain  asymptotic  conditions 
(e.g.  (6.14)  )  causes  some  complication.   But  the  integrals 
to  be  evaluated  in  a  collision  calculation  are  of  precisely 
the  same  kind  as  occur  in  the  binding  energy  calculation  of 
the  corresponding  bound  state.   This  indicates  that  scattering 
cross  sections  can  be  obtained  by  an  amount  of  numerical  work 
comparable  to  that  involved  in  bound  state  calculations. 

On  the  other  hand,  the  variation  iteration  method 
(  Part  II  ),  although  yielding  strikingly  good  results  in 
two  body  collision  problems,  appears  to  be  unsuited  for  the 
treatment  of  more  complicated  scattering  processes.   For  • 
we  have  seen  that  the  theoretically  most  satisfactory  procedure 


104 


(  section  II. 2.  )  is  completely  useless  from  a  practical 
point  of  view,  while  other  methods  (  section  11,3.  )  are 
still  extremely  complicated  and  of  doubtful  value.   Of  coursp 
one  cannot  rule  out  the  possibility  that  some  better  iteration 
procedure  exists.   However  any  integral  formulation  increases 
the  dimensionality  of  the  integrals  to  be  found.   Thus,  unless 
some  unforeseen  simplification  can  be  made,  we  do  not  believe 
that  any  such  method  can  be  as  useful  as  the  Rayleigh-Ritz 
me  thod • 

While  this  paper  falls  far  short  from  providing  a 
complete  solution  of  the  problem  of  composite  nuclear 
collisions,  it  is  hoped  that  it  has  made  some  contribution 
towards  paving  the  way  for  a  system;atic  treatment  of  such 
processes. 
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